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Abstract 


A new methodology for multiratc samplcd-data control system design based on (1) a new generalized 
control law structure, (2) two new parameter-oplimization-based control law synthesis methods, and (3) a 
new singular-value-based robustness analysis method is described. The control law structure can represent 
multirate samplcd-data control laws of arbitrary structure and dynamic order, with arbitrarily prescribed 
sampling rates for all sensors and update rates for all processor states and actuators. The two control law 
synthesis methods employ numerical optimization to determine values for the control law parameters to 
minimize a quadratic cost function, possibly subject to constraints on those parameters. The robustne- 
analysis method is based on the multivariable Nyquist criterion applied to the loop transfer function for the 
sampling period equal to the period of repetition of the system’s complete sampling/update schedule. The 
complete methodology is demonstrated by application to the design of a combination yaw damper and modal 
suppression system for a commercial aircraft. 
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I. Introduction 


The original objective for this project was to demonstrate a new algorithm for synthesizing multiraic samplcd- 
data control laws by application to a representative aircraft control problem. That algorithm, developed in 
connection with another research effort supervised by the Principal Investigator and based on a finite-time 
quadratic cost function, eventually proved unsuitable for the aircraft control problem. To complete this project 
we therefore developed a new multiratc control law synthesis algorithm, based on an infinite time quadratic 
cost function, along with a new method for analyzing the robustness of multiratc systems, and applied both to 
the aircraft control problem 

The following is a complete list of the contributions of this project: 

1 . A new generalized multiratc samplcd-data control law structure (GMCLS) was introduced. Features of 
this structure include an arbitrary dynamic order and structure for the processor dynamics; and sampling 
rates for all sensors, update rates for all processor states, and update rates for all actuators that can be 
selected independently, (discussed in Section II) 

2. A new infinite-time-based parameter optimization multirate sampled-data control law synthesis method 
and solution algorithm were developed, (discussed in Section HI) 

3. A new singular-value-bascd method for determining gain and phase margias for multirate systems was 
developed, (discussed in Section IV) 

4. The finile-time-based parameter optimization multiratc samplcd-data control law synthesis algorithm 
originally intended to be applied to the aircraft problem in this project, was instead demonstrated by 
application to a simpler problem involving the control of the tip position of a two-link robot arm. 
(discussed in Sections III and V) 

5- The GMCLS, the new infinite-time-based parameter optimization multirate control law synthesis method 
and solution algorithm, and the new singular-value based method for determining gain and phase 
margins were all demonstrated by application to the aircraft control problem originally proposed for this 
project (discussed in Section VI) 

These five contributions are discussed in order in the following sections of this report. The first three sections 
arc in a summary form only and the reader is referred, for details, to preprints of journal papers in the 
appendixes. The next two sections present applications of the parameter optimization techniques. The final 
two sections present our conclusions and suggest topics for future research. 
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II. Tiie Generalized Multirate Sampled-Daia Control Law Stricture 

A key point often ignored by the developers of multirate samplcd-daia control law synthesis methods is that, in 
order for any such method to be practically useful, it must provide the control law designer with the flexibility 
to independently choose the sampling rate for every sensor, the update rate for every processor state, and the 
update rate for every actuator. Such flexibility is frequently essential for efficient utilization of real-umc 
control hardware, and for systems that include distributed processing and/or utilize sensors that provide only 
discrete-time signals at fixed sampling rates [1], In this section we present a general-purpose, multirate 
sampled-data control law structure (GMCLS) that provides that flexibility. 

To understand the GMCLS, it is necessary to establish a certain notation regarding the scheduling of sampluig 
and update activities for a multirate system. Figure 1 shows an example of the time lines for the sampling and 
update activities of a multirate system. We define the shortest time period (STP) as the greatest common 
divisor of all of the sampling, update and delay periods; and we define the basic time period (BTP) as the least 
common multiple of all of the sampling, update and delay periods. We reserve the symbol T to represent the 
STP, and the symbol P to represent the (integer) number of STP’s per BTP. Finally, we frequently make use 
of a doubly-indexed independent (lime) variable, so that, for example, x(m.n) represents x at the start of the 
(/t+l)/A STP of the (m+\)th BTP, for m= 0,1, ... and «=0,1 P- 1. 

A block diagram of the GMCLS is shown in Figure 2. y represents the incoming, noise-free, continuous- 
time sensor signal; v is the discrete-time sensor noise signal; and u is the continuous-time control signal. The 
sampling period of the one sampler is the STP of the complete system’s sampling/updale schedule. The delay 
blocks are one-STP delays; and the ZOH block represents a zero-order hold. 


I'imc Lines for Sampling/Update Activities: 
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Figure 1 Example Multirate Sampling/Updatc Schedule 
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Figure 2 Generalized Multiratc Samplcd-Dala Conirol Law Structure 


A key feature of the GMCLS is its use of the switching matrices, S y (n), S 2 (n), and S u (n), for n=0,l P-l, 

to represent the variations in the sensor sampling, processor stale update, and control update activities, 
respectively. We define a switching matrix as a binary, diagonal matrix. S y (n) is the switching matrix that 
describes the sensor sampling activities at the start of the (n+\)lh STP (of every BTP). If the ilh diagonal 
element of S y (n) is 1 , then the ith sensor’s signal is sampled at the start of the {n+\)th STP of every BTP, and 
the sampled value, with the sensor noise v added, is immediately stored as the ilh element of y. If the ith 
diagonal element of S y (n) is 0, then the same element of y is simply held at those iastants. The update activities 
for the processor state vector z and for the actuator told state vector u, in Figure 1 , are similarly represented 
by the switching matrices S z (n), and S u (n), respectively, for n=0,l, ... ,/M. 

For a detailed discussion of the GMCLS see [ 1 ). The key points arc: 

1 . The switching matrices S y (n), S z (n), and S u (n) are completely determined by the system’s sampling and 
update activities schedule. 

2. The only unknowns are the processor matrices A z (n), B z (n), C 2 (n), and D z (n) 

3 . The dynamic order of the processor dynamics (i.e., the dimension of z) is arbitrary. 

For design purposes, the implications of these points are the following: 

1 . The GMCLS provides complete flexibility with regard to the selection of sampling rates for all sensors, 
update rates for all processor states, and update rates for all actuators. The single constraint is that the 
ratio of all sampling, update and delay rates must be rational, so that the complete sampling/updatc 
schedule is periodic. 

2. The GMCLS provides complete flexibility with regard to the dynamic order and structure of the control 
law; i.e., the input-output dynamics of virtually any multiratc samplcd-data control law of practical 
interest can be realized with the GMCLS. 
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3. Apart from the (significant) problem of choosing sampling and update rates, the GMCLS reduces the 
control law synthesis problem to one of determining the processor matrices A z (n), B z (n), C z (n), and 
D z (n), for n=0,l P-l. 

For the purpose of numerically determining A z (n), B z (n), C z (n), and D z {n) it is convenient to represent the 
GMCLS in the following stale model form (see [1 ] for details): 


c{m,n+\) = A c (n)c(m,n) + B c (n)y(m,n ) 
u(m,n) = C c (n)c(m,n) + D c (n)y(m,n) 

where 


c(m,n) = [z(m,n) y(m,n) u(m,n)] T 


A c (n) 


\l-S z (n))+S z (n)A z (n) 

0 

S u (nX^ z (n) 


S z (n)B z (n)[l-S y (n) \ 

l-S y (n) 

S u (n)D z (n)[l-S y (n)] 


0 

0 

l-S u (n) 


B c (n) 


~S z (n)B z (n)S y (n)~ 

S y (n) 

_S u (n)D z (n)S y (n)_ 


C c (n) = (S u (n)C z (/t) S u (n)D»[/-S y (n)] /-S„(n)] 
D c (n) = [S u (n) D z (n) S y (n)]y(mji) 
with «(t) = M(/n/i) for all l on [(mP + n)T, (mP + n + 1)7). 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 


The compensator parameters, A z (n), B z (n), C z (n), and D z (n), can be separated from the sampling schedule, 
S u (n), S y (n), S z (n), in an output-feedback representation of the GMCLS. Assuming a discretized model of 
the plant dynamics of the form 


p(m,n+ 1) = A p p(m,n) + B p u(m,n) + E p v(mji) (8) 

y(m,n) = C p p(m,n) + F p w(m,n) (9) 

where v and w represent represent process and measurement noise, respectively, we can rewrite the closed 
loop system in the output feedback form 
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( 10 ) 

( 11 ) 
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Now ihc compensator matrices can be factored as follows. 


DM) CM) 
L BM) AM) 


= SM) 


DM) CM) 

L BM) AM) J 


S2(n) + SM) 


(13) 


where Sj(/t), ^(n), and Sj(n) are functions of S u (n)< S y (n), and S z (n). 

Equation (13) is important because it allows us to separate the unknown compensator parameters A z (n), B z (n), 
C z (n), and D z {n) from the known sampling schedule. 


In the following section we will introduce two synthesis algorithms that can be used to determine the optimum 
compensator parameters A z (n ), B z (n), C z (n), and D z (n). 


III. Parameter Optimization Control Law Synthesis Methods 

There arc five well-recognized techniques for synthesizing multirate control laws: successive loop closures, 
pole placement, singular-pcrturbation-based methods, LQG Optimal methods, and parameter optimization 
methods. 

The advantages of successive loop closures are that it is easy to use, that it can (conceivably) be used to 
synthesize control laws of arbitrary dynamic order and structure, and that it is particularly effective in 
applications where the control loops are not strongly dynamically coupled. Its disadvantage is that its one- 
loop-at-a-time approach cannot fully account for all dynamic coupling between control loops. 

The problem with pole placement is determining where the closed-loop poles should be placed. It is a 
particularly difficult problem in the multirate (as compared to the single-rate) case because the STP-to-STP 
dynamics of multirate systems are periodically time-varying [2]. Only the BTP-to-BTP dynamics of multirate 
systems are time-invariant, and it is the poles of those dynamics that are assigned by pole-placement. In 
applications, determining desirable BTP-to-BTP closed-loop poles for a typical multirate system is difficult 
because the BTP of its sampling/updatc schedule will typically be longer than many of its desired closed-loop 
characteristic times. 

Singular-perturbation-bascd control law synthesis methods amount to successive loop closures prefaced with a 
coordinate transformation to separate the full control law synthesis problem into two or more dynamically 
decoupled control law synthesis problems of different time scales. A complete decoupling requires changes in 
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not just the state coordinates, but in the input and output coordinates as well. Such a decoupling is not 
possible in the multirate case because the input and output coordinates represent physical sensor and actuator 
signals destined to be samplcd/updated at different rates. 


The advantages of the LQG optimal control law synthesis methods arc that stabilizing control laws arc 
relatively easy to obtain and that the control laws for all control loops arc synthesized simultaneously, taking 
full advantage of all dynamic coupling between the control loops. The disadvantages are that the dynamic 
order and structure of the control law is fixed, that stability robustness objectives arc difficult to achieve, and 
that the resulting control laws are periodically time-varying |2)-[3]. 

We favor parameter optimization methods for control law synthesis for mulliratc systems because they offer 
the principal advantages of the successive loop closures and LQG optimal synthesis methods. These 
advantages are that control laws of arbitrary dynamic order and structure can be synthesized, and that control 
laws for all control loops can be synthesized simultaneously, taking full advantage of all dynamic coupling 
between control loops. The disadvantage of parameter optimization methods is that a numerical search is 
required to determine the control law parameters. 

In this section we present two parameter optimization methods for synthesizing multirate control laws. Both 
utilize the GMCLS discussed in Section II. The first is based on a finite-time quadratic cost function while the 
second is based on an infinite-time quadratic cost function. Both methods solve the multirate compensator 
synthesis problem by using a gradient-type numerical search to find a set of compensator parameters that 
minimize a quadratic cost function. 

The multirate optimization problem is as follows. 

Given: 

1 . The plant dynamics represented by 

p (/) = A p p(t) + B pu u(l) + B pv v(/) ( 14) 

y(i) = C p p(t) (15) 

Here p is the plant state vector, u is the control input vector, y is the noise-free measurement output 
vector, and v is the noise input vector. 

2. The complete sampling and update schedule for the compensator. This amounts to specifying S u (n), 

S y (n), and S z (n). for n=0,l P-l . 

3 . The order for the processor dynamics (the number of elements in z in (3)). 

4. The desired structure (e.g., a diagonal structure) for the processor matrices. A z (n), B z (n), C z (n), and 
D z (n), for n=0,l, ... , P-1. 

5. The number of distinct sets of processor matrices and when they are active. The optimization algorithms 
allow A z (n), B z (n), C z (n), and D z {n) to be periodically time varying. The designer can specify equality 
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relations among the compensator matrices. For example, if a time invariant compensator is desired then 
the designer can specify that A z (0) = ^(U = - = A z (P-l), and similarly for B z ,C Z and D z . 

6. The power spectral density V of the process noise v (in (8)) 

7. The covariance W(n), for «=0.1, ... , P- 1, of the sensor noise w (in Fig. 2). w ,s assumed to be a 
periodically stationary, gaussian, purely random sequence, with period equal to the BTP of the 
sampling/updatc schedule. 

8 . The time //and non-negative definite weighting matrices Q and R for the performance index 


{ V 

r - 

T 



r 

> 

r / 

Pit) 


Q 0 


p(o 

di ^ 

K J 

m. 


. 0 R _ 


_«(/)_ 

> 


where £ is the expected value operator. 


In the finite time optimization problem //must be a multiple of the BTP of the sampling/update schedule. 
In the infinite lime optimization problem // and J m f imte time = lim J(ir) 

lf~> OO 


Find: 

A set of processor matrices, A z (n), B z (n), C z (n), and D z (n ), for n=0.1 />- 1, such that the performance 

index 


J ss = lim J(tj) 

l j- oo 

is minimized. 

This optimization problem can be solved using either the finite-time cost function or the infinite-time cost 
function. 

Salution Method Using the Finite-Time Cost Function. 

The finite-time optimization algorithm was developed in connection with another research effort supervised by 
the Principal Investigator. This method synthesizes the multirate compensator that minimizes /(//) for a finite 
tf. A detailed discussion of this method can be found in [ 1 ]. A summary of the solution procedure follows. 

1 . Determine closed-form expressions for the performance index and for its gradients with respect to 

the elements of the processor matrices A z (n\ B z (n), C z (n), and D z (n ) . for n= 0.1 P-l . 

2. Use a gradient-type numerical optimization algorithm to determine a set of processor matrices, A z (n), 
B z {n), C z (n), and D z (n ) , for n=0,l, ... ,P-1, that minimizes 
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3 . Obtain a steady-stale solution by re-optimizing for larger and larger until tj- gets to be large compared to 
all of the closed-loop system’s characteristic times. 

The advantage of this method is that with tj- finite, the cost function J(tj-) remains finite even if the compensator 
is destabilizing. The designer does not need to find a stabilizing compensator to start the optimization process 
as long as /y is small enough that 7(/y) does not exceed the numerical limits of the computer performing the 
optimization. 

The disadvantage of this method is that the closed- form expressioas that have been developed thus far for the 
performance index J(tj) and for its gradients with respect to the elements of the processor matrices arc very 
complex and computationally intensive. In addition, we encountered difficulties when applying this method to 
the aircraft control problem because our solution algorithm lacked provisions for automatic scaling of the 
control law parameters (i.c., the independent variables) during the numerical search. The sheer complexity of 
the finite-time performance index and gradient expressions prevented us from adding the automatic scaling 
provisions that would have allowed us to apply this method to the aircraft control problem. 

Solution Method Using Infinite -Time Cost Function 

Instead of modifying our existing finite -time-based algorithm to alleviate the scaling problem discussed in the 
previous paragraph, we chose to develop a new infinite-time-based multirate sampled-data control law 
synthesis method, based on corresponding developments for single-rate systems by Mukhopadhyay [4], for 
which much simpler performance index and gradient expressions are easy to derive. For a complete 
description of that method, and the solution algorithm we developed to implement it see [5]. A summary of 
the solution procedure follows. 

1. Find an initial stabilizing guess for the processor matrices A z (n) y B z (n),C z (n ), and D z (rt) % for 

n= 0,1 P-[. The finite -time solution algorithm requires an initial stabilizing compensator because 

J ss is infinite when the closed loop system is unstable. From our experience, many multirale problems 
can be stabilized using successive loop closures. The aircraft problem was open loop stable, and so 
determining a stabilizing compensator was trivial. 

2. Determine the necessary conditions (given in [5]) for the processor matrices, A z (n ), B z (n)> C 2 (n) y and 
D z (n), for n=0,l, ... ,/M to minimize J ss . These are represented by three sets of coupled matrix 
equations. Two sets are Lyapunov equations, one governs the steady state covariance of the plant and 
control states, and the other governs a Lagrange multiplier. The third represents the gradient of J ss with 
respect to the compensator parameters. 

3. Use a gradient-type numerical search to solve the necessary conditions and determine a set of processor 
matrices, A z (n) y B z (n) y C z (n) y and D z (n) f for n=0,l, ... ,/M, that minimizes J ss . 

The advantage of this method is that the gradient of J ss with respect to the compensator parameters is easy to 
evaluate via the necessary conditions. For a given problem, the infinite-time optimization algorithm typically 
requires fewer computations to find the optimum compensator parameters than does the finite-time 
optimization algorithm even when both algorithms are initialized with the same stabilizing compensator. 
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Even .hough U* finite-time and mfiniie-iime based soludon algonduns can dc.erm.ne opdmum compcnsaiur 
parameters, toe is no guaran.ee to Ihc des.gn will be robust. In .he follow, ng section we presen, a method 
Tor analyzing the robustness of a multirale control system. 


IV. Gain and Phase Margins for Multirate Systems Using Singular-Vali es 

There are many established methods for synthesizing multiratc compensators, see Section ffl. but surprisingly 
few methods for analyzing the robustness of these systems. Current robustness analysis methods rely 
pnne. pally on the transfer function of the system. A multiratc transfer function, in the traditional sense, does 
not exist, because mult.ratc systems arc periodically time varying. Without modification, established single 
rate analysis methods cannot be applied directly to multiratc systems. 

As par, of Ihis project, we developed an approach for extending Ac nyquist eri.enon and singula, value 
analyse ,o mult, rate and periodically time vary.ng syslems. For a detailed dtscusston of to approach 
including appltcauon of structured singular value robustness analysts to multi, ale systems, see [6j. In this 

secuon we present a summary of the important ideas from to paper used to calculate gain and phase margins 
of multirate systems using singular values. 

As we saw in Section II. a multirate compensator can be modeled as a I, near periodically time varying system 
U)-(2). Equations (1 M2) from Section II can be written as 


1 ) - A c (n)c(m,n) + B c (n)y(m,n) ^ j 7 j 

u{m y n) = C c (n)c(msi) + D c (n)y(m.n) ( l 8 ) 

This system ( 1 7)-( 18) can then be transformed to an equivalent single-rate system (ESRS) by repealed 
application of (17)-(18) over the BTP |7J. The ESRS has the forni: 


c(m+l ,0) = A e c(m,0) + B e y(m, 0) 
u(m,0) = CeC(m, 0) + Djy(m,0) 


where y(m, 0) = 

” y(m, 0) ' 
>(w,l) 

u(m, 0) = 

” u(m, 0) 
u(n i.l) 


1). 




(19) 

( 20 ) 

( 21 ) 


The transfer function for the ESRS is 


( 22 ) 


u(z P ) = Gp(z p )y(z p ) 
where Gp(z p ) = C e (Iz p - A e ) ] B e + D e 


( 23 ) 



For a detailed discussion of the ESRS, see [6], The key points arc: 

1 . The ESRS is a time invariant single-rate system with a sampling period of one BTP and the unique 
property that the inputs arc time correlated and the outputs arc time correlated. 

2. In general G P (z p ) has a very complicated form, but it can be shown that if the system is time invariant 
with G(z) equal to a constant, then G P (z p ) will also be constant and block diagonal with G(z) on die 
diagonal. 

3. The ESRS allows us to manipulate time invariant and periodically time varying systems (e.g. multirate ) 
as if they were both time invariant. The state space or transfer functions descriptions can be used to 
calculate input-output relations for systems in series or in a feedback loop just as in classical control [8], 
For example, to calculate the ESRS of a multirale compensator in scries with a time invariant plant, wc 
would calculate the ESRS of the plant and compensator individually and then combine them using block 
diagram arithmetic. 

4. Kono [9] has shown that if the ESRS is stable then the mullirate system from which it was derived will 
be stable. 

5. Single-rate robustness analysis techniques can be applied to the ESRS as long as the results are 
interpreted in light of the fact that some of its inputs and outputs are time correlated. 

Generalized gain and phase margins for the ESRS (and equivalently the mullirate system) can be calculated 
using singular value analysis. If wc assume a plant uncertainty of the form 

Actual ~ G(z)^ <)mlna ikt^ (24) 

then the ESRS plant uncertainty has the form 

C, P^ P '> Actual - G p( zP )Nommal( lce,Q )p ( 25 ) 

(kefi) P = diagt^ 0 , ke > Q , .... ke f®] with P blocks 

[Recall that if H(z) is constant, H P (z p ) is block diagonal with H(z) on the diagonal.] 

The multirate system is guaranteed to remain stable whenever 

Q(l + G P (z p )) on the nyquist contour (26) 

Traditional gain margins can be obtained by setting 0 = 0 and solving (26) for k. Phase margins can be found 
by selling k = 0 and solving (26) for 0. 

As with most singular value robustness analysis methods, the it and 0 found using (26) are conservative. If, 
however, ke> Q is diagonal, the conscrvativcncss associated with (26) can be reduced by diagonally scaling 
Gf>(z p ). Wc used Osborne’s method of preconditioning matrices to increase the lower bound for G P (z p ) and 
thus to improve our estimate of the gain and phase margins. 
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V. Application of the Finite-Time-Based Parameter Optimization algorithm to 

a Two Link Robot arm Control Problem 

The original proposal for this project called for the finite-time-based parameter optimization muluraie sampled- 
data control law synthesis method of Section HI to be applied to an aircraft control system des.gn problem. 
That method and a solution algorithm to implement it had been previously developed as part of another 
research effort supervised by the Principal Investigator. Due to the solution algorithm's lack of adequate 
provisions for automatic scaling of the control law parameters (i.e., the independent variables) during the 
numerical search, we were not able to apply it successfully to the aircraft control problem. We maintain, 
however, that the problems we encountered with it were a consequence of problems with the solution 
algorithm and are not necessarily indicative of problems with the synthesis method. 

In this section we therefore present an application of the fimte-time-bascd mulurate samplcd-data control law 
synthesis method to a two-link robot arm (TLA) control problem. The robot arm application demonstrates the 

utility of the method without being so poorly conditioned that automatic scaling of the control law parameters 
was required during the numerical search. 

The two-link robot arm system we dealt with is shown in Figure 3. The first link is long and massive, for 
large-scale slewing motions. The second is short and lightweight so high-bandwidth control of the tip position 
can be achieved with a relatively small motor at the second joint. The pin joint, rotational spring, and 
rotational damper at the midpoint of the first link model flexibility in that link. The control inputs are the motor 
torques, T, and T 2 . The measured outputs are the joint angle 6 and the tip position 5. The spring constant (k) 

and damping coefficient (b) values (in Fig. 3) yield an open-loop vibration mode with a 10 Hz natural 
frequency and 1 % damping. 



L x 0.5 kg 0.5 m 

L 2 0.5 kg 0.5 m * = 37.33 N/rad 

L 3 0.04 kg 0.2 m b = 0.012 N- s/m 

The natural frequency of the vibration mode is 10 hz. 
Injxits: Torques T\ and T 2 
Outputs: 0 and 8 

Figure 3 Two-Link Robot Arm System 
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We used the finite-time-based multirate control law synthesis method of Section III to synthesize multiratc 
sampled-data control laws for this system. Our performance objective was high-bandwidth control of the tip 
position 8, and it is intuitively clear that this can best be accomplished, given a fixed real-time computation 
capability, by trading low-bandwidth control at T\ for high-bandwidth control at 7 2 . Thus, for an 8-to-l 
control bandwidth ratio, we chose the sampling/update rate for 5 and 7 2 to be 8-limes faster than that for 0 and 
T\. For comparison purposes, we also designed corresponding analog and single-rate sample-data control 
laws. 

TLA CquuqI La w s 

For the TLA system, the tip position (8) responses to a commanded step change in the tip position obtained 
with the analog, single-rate and multirate control laws we synthesized are shown in Figure 4. Sec ( 1 J for 
additional results and details. A summary description of those designs follows. 

LQR Analog Design The LQR Analog response was obtained with an analog LQR (full state feedback) 

control law that is optimal with respect to a quadratic performance index that yields 0.7071 damping 
(£j = £ 2 = 0.7071) and an 8-to-l ratio of characteristic frequencies (co n2 /w nl = 8) for the two closed-loop 

modes. 

Third-Order Analog Successive Loop Closures Design The Third-Order Analog Successive Loop Closures 
response was obtained with a successive loop closures control law that consisted of a single lead network from 
0 to Tj, and two identical cascaded lead networks from 8 to T 2 . The gains, and zero and pole locations were 
chosen to yield dominant closed-loop poles coincident with those obtained with the LQR Analog control law. 

Third-Order Multirate Tustin Design The Third-Order Multirate Tustin response was obtained with a 
multirate sampled-data control law obtained by discretizing the lead compensators of the Third-Order Analog 
Successive Loop Closures design using Tustin’s method [10]. The O-to-Tj control-loop sampling/update rate 
is a factor of 8 times the characteristic frequency of the lower-frequency closed-loop mode from the Third- 
Order Analog Successive Loop Closures design; and the 8-to-T 2 sampling/update rate is the same multiple of 
the characteristic frequency of the higher-frequency closed-loop mode from the Third-Order Analog 
Successive Loop Closures design. 

Optimized Third-Order Multirate Tustin Design The Optimized Third-Order Multirate Tustin response was 
obtained with a control law synthesized by the finite-time-based multirate sampled-data control law synthesis 
method of Section III. This control law is the Third-Order Mullirale Tustin control law, but with its gains and 
its pole and zero locations optimized to minimize the same performance index as is minimized by the LQR 
Analog control law. 

Analog Third-Order Design The Analog Third-Order response was obtained with a third-order, generally- 
structured, analog control law synthesized using Ly’s Sandy algorithm [ 1 1 ]-[ 1 2] to minimize the same 
performance index as is minimized by the LQR Analog control law. 
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Multiraie First-Order, Second-Order & Third-Order Designs The Mullirate First-Order, Second-Order, and 
Third-Order responses were obtained with mullirate, generally-structured, sampled-dala control laws 
synthesized by the finile-time-based multirate samplcd-data control law synthesis method of Section HI to 
minimize the same performance index as is minimized by the LQR Analog control law. The sensor sampling 
and actuator update rates are the same as in the Third-Order Mullirate Tustin control law. In the First-Order 
case, the update rate for the one processor stale is the same as the faster sensor-sampling/actualor-updatc rate. 
In the Second-Order case, one processor state is updated at the faster rate and the other at the slower rate. In 
the Third-Order case, two processor states are updated at the faster rate and one is updated at the slower rate. 

Single-Rate Third-Order Design Finally, the Single-Rate Third-Order control law response was obtained 
with a single-rate, generally-structured, sampled-data control law synthesized by the finite-time-based multirate 
sampled-data control law synthesis method to minimize the same performance index as is minimized by the 
LQR Analog control law. Its single sampling/update rate was chosen to require the same average number of 
computations per unit lime for real-time operation as is required for real-time operation of the Multirate Third- 
Order control law. 

Conclusions 

The TLA results in Figure 4 demonstrate some of the benefits of multirate control. For example, the tip 
position overshoot (6) with the multirale compensator is much less than with its equivalent single-rate 
counterpart. But more importantly, the results demonstrate that the finile-time-based multiraie sampled-data 
control law synthesis method can be used to synthesize mullirate control laws of arbitrary structure and 
dynamic order, with arbitrarily selected sampling rates for all sensors, and update rates for all processors states 
and actuators. The third-order mullirate compensator, for example, uses two different update rates for the 
processor states, inputs and outputs, and a general compensator structure with full coupling between inputs, 
outputs and processor states of different rates. 
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VI. APPLICATION OF THE INFINITE-TIME-BASED PARAMETER OPTIMIZATION ALGORITHM 

to a Yaw Damper and Modal Suppression System for a Commercial aircraft 

A practical application of mulliratc control can be found in aircraft. The limited computational resources of 
aircraft dictate that their control systems must function efficiently. Mulliratc control allows the designer to 
efficiently allocate these resources by trading slow sampling and update rates in control loops associated with 
low-bandwidth control functions for fast sampling and update rates in control loops associated with high- 
bandwidth control functions. In this section we consider a particular application of mullirate control: a 
combination yaw-damper and modal suppression system for a commercial aircraft. 

In the interest of weight reduction for fuel efficiently, aircraft are being constructed with less structural rigidity. 
Structural vibration modes can be excited in such aircraft by wind gusts or by movements of control surfaces. 
These vibrations affect not only the structural integrity of the fuselage but also passenger ride quality. In the 
lateral direction, such vibrations are often induced by rudder activity associated with the yaw-damper. A 
“modal suppression system” can be added to the yaw-damper loop to suppress these vibrations. The modal 
suppression system would traditionally be designed by successive loop closures. 

In this section we describe the design of a multirate combination yaw-damper and modal suppression system 
for a commercial aircraft using the infinite-lime-based multirate compensator synthesis algorithm and 
robustness analysis technique discussed in Sections III and IV. For comparison purposes we also designed 
corresponding analog and single -rate sample-data systems. 

The goal for each compensator design was to increase the damping of the dutch-roll mode to 0.6, and to 
decrease the covariance of lateral accelerations at the nose and aft of the airplane, particularly those components 
associated with low frequency flexible modes. The performances of the compensators were compared by 
comparing the closed loop dutch-roll damping, the covariances of lateral accelerations at the nose and aft of the 
aircraft due to a unit covariance gaussian white noise disturbance, and the PSD plots of the lateral accelerations 
at the nose and aft of the aircraft for either a white noise disturbance (analog designs) or a gust pulse 
disturbance (samplcd-data designs). 

Open Loop Aircraft 


A block diagram of the airplane model is shown in Figure 5. The lateral dynamics model consists of 4 rigid 
body modes (heading, spiral, dutch roll and roll) and 1 1 flexible modes. Actuator/power control units for the 
aileron and rudder are modeled as second-order lags. 


(*($) Rudder ~ + 3 ^)( s + 35) 


rf v 20(25) 

^ViAileron ~ ( s + 20)(r + 25) 


(27) 


The lateral gust disturbances are filtered by a second-order Dryden gust model 


G(s) = 


17.496* + 2.1617 
21.836* 2 + 9.3458* + 1 


(28) 
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Figure 5 Block Diagram of the Open Loop Airplane Model 


The poles associated with the spiral and heading modes were compensated with static gain feedback before the 
yaw-damper/modal suppression systems were designed, because these modes, which laid close to the origin 
and were controllable with the rudder, created numerical difficulties for Sandy (the optimization program used 
to design the Fourth-Order Analog compensator discussed in later in this section). The spiral mode was 
compensated by feeding back roll and roll-rate to the aileron. Heading was compensated with heading to 
rudder feedback. In what follows we refer to the airplane model with spiral and heading modes compensated 
as the uncompensated airplane (no dutch-roll compensation). 

The lateral accelerations of the uncompensated airplane are measured by Nynose and Nyaft. The PSD plots of 
lateral accelerations for the uncompensated airplane are shown in Figure 6. A yaw-damper/modal suppression 
system should reduce the total area under this curve (covariance of lateral acceleration). In particular, it should 
reduce the peak at = 0.5 Hz (near the dutch-roll mode) and the peaks between 3 Hz and 6 Hz (low frequency 
flexible modes). Values of the dutch-roll damping, and the Nynose and Nyaft covariances for the 
uncompensated airplane are given in Table 1 . 
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Figure 6 PSD of Nynose and Nyafl for Uncompensated Airplane with Unit 
Covariance Gaussian White Noise Lateral Disturbance 


Table I Results for Analog Designs with a Unit Covariance Gaussian White Noise 
Lateral Disturbance 


Design 

Dutch-Roll Damping 

Nynose Cov. 
(ft 2 /scc 3 ) 

Nyafl Cov. 
(ft 2 /scc 3 ) 

Uncompensated 

0.08 

5.1 

21.8 

Analog Yaw-Damper Only 

0.6 

5.0 

6.1 

LQR Analog 

0.6 

2.4 

3.1 

Fourth-Order Analog 

0.55 

2.5 

2.4 


Analog Yaw-Damper/ Modal Suppression System Designs 

Three analog compensators were designed: a yaw-damper only system, a full state feedback yaw- 
damper/modal suppression system, and a fourth-order yaw-damper/modal suppression system. PSD plots of 
Nynose and Nyafl for the analog designs are shown in Figure 8; Nynose and Nyafl covariances are 
summarized in Table 1. These analog designs provide a base line for comparison with the sampled-data 
designs and were used to determine appropriate values for cost function weighting matrices. Following is a 
summary of these designs. 
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Analog Yaw-Damper Only Design The yaw-damper only design uses sialic feedback from y to 8 r , using a 
gain k-damper- We chose k damper such that the dutch-roll damping was 0.6 using classical root locus. While the 
peak on the PSD plot associated with the dutch-roll mode (=0.5 Hz) has been reduced significantly from the 
uncompensated case, the peak near 3 Hz has increased (see Fig. 8). This is the problem with using static gain 
feedback. As you "press” on one peak of the PSD another “pops” up due to the input coupling between the 
dutch-roll and low frequency flexible modes. 


LQR Analog Design The LQR design uses full state feedback to improve the dutch-roll damping and reduce 
the covariance of Nynose and Nyafl. The compensator was designed to minimize the following cost function. 


Jss = lim E 


j Nynose 
_ Nyaft 


T 


0.001 0 
0 0.004. 


Nynose 
. Nyaft . 


+ l .65 


(29) 


Weighting matrices for (29) were chosen such that the covariances of Nynose and Nyaft were reduced from 
the yaw-damper only case by the same percent, and the dutch-roll mode had a damping of 0.6. Figure 8 
shows that the LQR design significantly reduces the dutch-roll peak as well as the peaks associated with the 
flexible modes. 

Fourth-Order Analog Design The Fourth-Order Analog compensator is a yaw damper/modal suppression 
system designed using Sandy [11]. A block diagram of this compensator is shown in Figure 7. This design 
minimizes the same cost function as the LQR design (29) with the weighting on 8 f adjusted to achieve close to 

the desired 0.6 dutch-roll damping. An unexpected result is that the covariance at Nyaft for the Fourth-Order 
Analog design is actually better than for the LQR Analog design. This is a consequence of adjusting the cost 
function weighting matrices to achieve the desired the dutch-roll damping. 


Lateral Disturbance 


Uncompensated Airplane 


¥ 

Nynose 

Nyaft 


Fourth- Order Analog 
Compensator 


Figure 7 Block Diagram of Airplane with Fourth-Order Analog Compensator 
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Figure 8b PSD of Nyafl for Airplane with Analog Compensators for Unit 
Covariance Gaussian White Noise Lateral Disturbance 
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Sampled-Data Yaw-Dampcr/Modal Suppression System Designs 


Three sampled-data compensators were designed: a single-rale yaw-damper only system, a fourth-order 
multirate compensator and a fourth-order single-rate compensator. Both fourth-order compensators were 
synthesized using our inTmite-time-based multirate control law synthesis algorithm to minimize the same cost 
function as the L QR Analog design. 

The sampled-data compensator designs were based on a maximum sample/update rate of 50 Hz. This is 10 
times the rudder actuator roll off frequency and 8 times faster than the fastest flexible mode which contributes 
significantly to the PSD of the lateral acceleration. This sample rate is close to the slowest practical sample rate 
which could be used. 

PSD plots for the sampled-data designs were generated using a gust pulse (a rectangular pulse) at the 
disturbance input, as opposed to the gaussian white noise used for the analog designs. For the analog 
designs, the PSD plots were based on transfer functions from the disturbance input to Nynose and Nyaft . 
Multirate compensators are periodically time varying so that transfer functions for them, in the traditional 
sense, do not exist. For this reason, we used the gust pulse disturbance input to generate the PSD plots for the 
sampled-data designs. 

The gust pulse input PSD has a connection to the white noise input PSD. For a time invariant continuous 
system, the PSD plots generated using either gaussian white noise or a continuous impulse input are exactly 
the same. This is because the Fourier transform of the impulse respoasc is the same as the bode plot, and the 
PSD of gaussian white noise is a constant. If a continuous system, given by Jc(/) = Ax(i) + Bu{i) y is such that 

T p 

B = u Je AT Bd x (30) 

where u can be selected arbitrarily and T p is much shorter than the observation time, then a continuous impulse 
can be approximated by a pulse of duration T p and magnitude u. For the airplane problem addressed in this 
project, (30) is satisfied for 

T p = 0.02 seconds and u = 50 ft/sec. (3 1 ) 

PSD plots of Nynose and Nyaft for the sampled-data designs are shown in Figure 1 1 . Nynose and Nyaft 
covariances for these designs are summarized in Table 2. Following is a summary of the sampled-data 
designs. 

Single-Rate Yaw-Damper Only Design The Single-Rate Yaw-Damper Only design is similar to the Analog 
Yaw-Damper design except that a sampler is used at the output y and a zero order hold is used at the input b r 
Both the sampler and zero order hold operate at 50 Hz. The performance of the sampled-data yaw damper is 
very close to that of the analog damper (Figs. 8 and 1 1). 
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Multirate Fourth-Order Design The muhiraie compensator is shown in Figure 9. It was designed to 
minimize the same cost function as the LQR Analog design with the weighting on 5 r adjusted to achieve the 

desired dutch-roll damping. The compensator uses two samphng/updatc rates. The rudder is updated and the 

lateral accelerations arc sampled at 50 Hz; y is sampled at a slower rate, 12.5 Hz, because it is composed 
primarily of the slow dutch-roll mode. 

Two of the processor states for this multirate compensator are updated at the fast rate, 50 Hz, and two are 
updated at the slow rate, 12.5 Hz. Initially a compensator was designed in which all of the processor states 
were updated at 50 Hz, but we found that there was no noticeable performance degradation if two of the 
processor states were updated at the slower rate. Slowing the update rate of these states reduces the number of 
computations required per unit time for real time implementation of the multirate compensator. 


Table 2 Results for Sampled-Data Designs with a 
Noise Lateral Disturbance 

1 . 

Unit Covariance Gaussian White 

Design 

Dutch-Roll Damping 

Nynose Cov. 

i A lyaft Cov. 



(ft 2 /scc 3 ) 

(ft 2 /scc 3 ) 

Uncompensated 

0.08 

5.1 

[ — 

21 X 

Single-Rate Yaw-Damper Only 

0.6 

4.3 

5.4 

Multirate Fourth-Order 

0.6 

3.6 j 

4.7 

Single-Rate Fourth-Order 

0.6 

3.5 ! 

4.7 



Figure 9 Block Diagram of Airplane with Multirate Fourth-Order Compensator 
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Single-Rate Fourth-Order Design The Single-Rate Fourth-Order compensator is shown in Figure 10. The 
sampling rate for this compensator is 28.6 Hz. That rate was chosen such that the number of multiplications 
required per unit lime for its real time operation is the same for the multirate compensators. The cost function 
used to design the single-rate compensator was the same as was used to design the Multirate Fourth-Order 
compensator. 



Figure 10 Block Diagram of Airplane with Single-Rate Fourth-Order Compensator 
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Single -Rate Yaw-Damper Only 
Multiratc Fourth-Order 
Single -Rate Fourth-Order 


Figure 11a PSD of Nynose for Airplane with Sampled-Data Compensators for 
Gust Pulse Lateral Disturt>ance 


Single-Rate Yaw-Damper Only 
Multirate Fourth-Order 
Single-Rate Fourth-Order 


Figure 1 lb PSD of Nyaft for Airplane with Sampled-Data Compensators for Gust 
Pulse Lateral Disturbance 





Phase 6 in Deg 


Gain and Phase Margins for Sampled-Data Design 

Gain and phase margins at the control input (6 r ) were evaluated for the Multirate Fourth-Order and Single-Rate 
Fourth-Order compensators using the robustness analysis methods of Section IV. Table 3 summarizes the 
traditional gain and phase margins for the these compensators. Figure 12 shows the region of guaranteed 
stability for simultaneous changes in k and 0 for both compensators. 


Table 3 Traditional Gain and Phase Margins 


Design 

Gain Margin (db) [0 = OJ 

Phase Margin (Deg) [* = Odb] 

Multirate Fourth-Order 
Single-Rate Fourth-Order 

[-3.8, 7.1] 
[-3.3, 5.5] 

±32* 

±2T 



Figure 12 Stability Region for Simultaneous Gain and Phase Uncertainly for Fourth- 
Order Sampled-Data Compensators 
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F.gures 8 and 1 1 show that the yaw-damper/modal suppression systems significantly decrease the covariance 
of the lateral acceleration at the nose and aft of the airplane while attaining the desired 0.6 dutch-roll damping. 
It should be no surprise that the analog compensators out performed the sampled-data compensators because 
the sampled-data compensators were designed using a slow sampling rate. Still, both fourth-order sampled- 
data compensators reduce the peak accelerations of Nynose and Nyafl by 1 75% and 50% respectively over the 
yaw-damper only systems. The performance of the Single -Rate Fourth-Order compensator is nearly as good 

as that of the multirate compensator, but, for input gain and phase uncertainty, the multirate compensator is 
more robust than the single-rate compensator. 


VII. Summary and Conclusions 

In this report we have presented a methodology for designing mullirate control systems. have introduced 
the Generalized Multirate Control Law Structure (GMCLS) which allows complete flexibility with regard to 
die dynamic order and structure of the control law, and with regard to die sampling rales for all sensors and the 
update rates for all processor states and actuators. We have presented two parameter optimization multirate 
control law synthesis algorithms, one based on an infinite-time cost function and the other based on a finitc- 
umc cost function, which can be used to find optimum values for the GMCLS parameters. We have presented 
a technique for determining gain and phase margins for multirate systems. Finally, we have demonstrated our 
methodology by applying ,t to the design of a two link robot arm control system and to the design of a 
combination yawnlamper and modal suppression system for a commercial aircraft. The application to the 
aircraft control problem, in particular, demonstrates that the methodology can be applied to design problems of 
a scale that one might expect to encounter in practice. 


VIII. Suggestions for Future Research 

The results presented here demonstrate a methodology for multirate digital control system design that is 

applicable to practical problems. Before this methodology can be routinely appl.ed in practice, however, the 
following need to be developed: 

1 . A means for directly synthesizing robust mullirate control laws. 

2. Numerical optimization algorithms incorporating auto-scaling of the independent variables and other 
features that more effectively deal with the practical difficulties of parameter optimization applied to 
multirate control law synthesis. 

With regard to direct synthesis for robustness, there are several possibilities. One would add the multiple- 
plant-condition design for robustness ideas of Ly [ 10]-[ 1 1 ]. A second would add direct nonlinear robustness 
constraints on the control law parameters during the numerical optimization. The latter approach has been 
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successfully applied by Mukhopadhyay [4] lo synthesize robust single-rate control laws by parameter 
optimization. 


In addition to theoretical work, a second major research effort in muluraic control needs to be directed toward 
experimental research. Now that a bonafide multirate control system design methodology has been developed, 
we strongly believe that further substantive progress in the field can best be made in conjunction with bonafide 
hardware applications of that methodology in the laboratory. 
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I INTRODUCTION 


Even in this age of fast, low-cost microprocessors there remain several important motivations 
for mul, irate sampling in sampled-data control systems. One ncx-d only consider large space 
structure control problems to realize lhal the cost, bulk, and weight of real-time computing 
hardware continues to be an important control system design issue. Multirate sampling provides 
the opportunity to allocate sampling rates, and thus real-time computing power, more efficiently. 
In two-time-scale control problems, for example, multirate sampling allows slow sampling in 
control loops associated with low-bandwidth control functions to be traded for fast sampling m 
those associated with high-bandwidth control functions. 


As with microprocessors, the costs of analog-, oaiigital and digital-, o-a„alog converters are also 

computalion-rate dependent. Mullirate sampling thus provides another opportunity to reduce 

hardware costs because the computation rales required of analog-lo-digital and digital-to-analog 

converters frequently depend upon their sampling rates. Multirate sampling can even be used to 

reduce the total number analog-, o-digital and/or digi,al-,o-analpg converters required by a system. 

by sample-dependent scheduling of multiple conversion tasks to a lesser number of conversion 
devices. 


A third "motivation" for multirate sampling is becoming increasingly important: sometimes 
multirate sampling is the only choice. This situation can arise when an apriori decision has been 
made to include in a system a sensor that provides a discrete-time signal at a fixed sampling rate. A 
head position control system for a computer disk drive is a good example of such a system. The 
disk head, which is suspended atop the rotating disk, includes a sensor that reads the head position 
directly from certain diametrically-spaced segments on the disk. The sensor's sampling rate is thus 
fixed by the disk's rotation speed. To increase the control bandwidth beyond that dictated by that 
sampling rate, a second, faster-rate sensor must be added. 


A key point often ignored by developers of multirate control law synthesis methods is that 
these motivations for multirate sampling dictate also certain flexibilities required to meet the needs 
of engineering practice Specifically, multirate control law synlhesis methods, to meet Ihe needs of 
engineering practice, must allow the sampling rates for all sensors, Ihe update rates for all processor 
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states, and the update rates for all actuators to be specified independently. The one generally 
accepted restriction, with regard to these rates, is that the ratio of all combinations of sampling and 
update rates must be rational, so that the complete sampling/update schedule will always be 
periodic. (We assume that all sampling and update events are synchronized to the same clock. 

The asynchronous case is treated elsewhere [1].) 

Time lines representing such a periodic sampling schedule are shown in Fig. 1. We define the 
Basic Time Period (BTP) of such a schedule as the least common multiple of all of its sampling and 
update periods. The BTP is the period of repetition of the sampling/update schedule. We define 
the Shortest Time Period (STP) as the greatest common divisor of all of its sampling and update 
periods. We reserve the symbol P to represent the (integer) number of STP’s per BTP, and we shall 
frequently use a double-indexing scheme for the independent variable so that, for example, x(m,n) 
represents x at start of the (n+l)th STP of the (m+l)th BTP, for m = 0,1 and « = 0 P - 1. 

There are five well-recognized methods for synthesizing multirate sampled-data control laws: 
successive loop closures, pole placement, the singular perturbation method, the LQG (linear 
quadratic Gaussian) method, and parameter optimization methods. Successive loop closures [21 is 
arguably the most important because it is the single one of the five that is widely used in industry. 
The advantages of successive loop closures are that its one-loop-at-a-time approach requires no 
new multirate synthesis techniques, and that the sampling/ update rate for each control loop can be 
specified independently. The problem with successive loop closures is that its one-loop-at-a-time 
approach cannot fully account for all dynamic coupling between control loops. 

Pole placement [3,4,5,6,71 for multirate systems has received considerable recent attention in the 
wake of reports on the capacity for periodically time-varying output feedback controllers to place 
closed-loop poles. In Ref. 3, for example, it is shown that given any controllable and observable 
continuous-time plant with m inputs, it is always possible to construct a periodically time-varying, 
purc-gain, output feedback control law that places the closed-loop poles arbitrarily, provided that 
the outputs are all sampled at a suitably chosen single sampling rate 1/To/ anc * that the inputs are 
updated at the rates N,/T 0 , N m /T c , where the N, are certain positive integers. 
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The problem with pole placement for multirate systems is the same as with pole placement for 
single-rate systems: how to determine where the closed-loop poles should be placed? It is a 
particularly difficult problem in the multirate case because multirate systems are periodically time 
varying [2,8]. The periodicity of multirate systems implies that their eigenstructure can only be 
defined based on their (time-invariant) BTP-to-BTP dynamics. Determining desirable closed-loop 
poles for a multirate system is typically difficult because the BTP of a multirate system is typically 
much longer than the characteristic times of many of its faster dynamics. 

Singular perturbation control law synthesis methods [9,10,11,12,13,14,15] were first developed 
for continuous-time control systems to take advantage of the multiple-time-scale dynamics that 
often occur in control systems. It would seem that an extension to multirate sampled-data systems 
should follow naturally, given that a principal motivation for multirate sampling has always been 
to take advantage of those same multiple time scales, but that has not been the case in practice. 

The problem is the singular perturbation method's inherant dependence on a coordinate 
transformation to separate the full control law synthesis problem into two (or more) dynamically 
decoupled control law synthesis problems of different time scales. Such a coordinate 
transformation is the first step in control law synthesis by the singular perturbation method. The 
state coordinates are easily decoupled because they represent only the plant's internal dynamics. 

The input and ouput coordinates cannot be so manipulated because they represent the plant's 
external sensor and actuator signals. Consequently, during the second control law synthesis step, 
when the control laws for the different-time-scale state vector components are synthesized 
separately, every control input vector element and every sensor ouput vector element remains 
coupled to every state coordinate so that, just as with successive loop closures, all dynamic coupling 
between control loops cannot be accounted for. 

Various schemes have been developed to circumvent this difficulty. None have been 
completely successful. In Ref. 13, for example, a state feedback control law is synthesized by the 
singular perturbation method, and the lack of a completely decoupling transformation gives rise to 
a requirement for the slow component of the plant state vector to be estimated between slow- 
sampler updates, and a requirement for every control input to be updated at every 
sampling/update instant. 
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The advantage of the LQG method 12,16,17,18] for multirate sampled-data control law synthesis 
is that the control laws for all control loops are synthesized simultaneously, taking into account all 
dynamic coupling between control loops. The disadvantages are the same as with the LQG method 
for continuous-time control law synthesis: that practical performance and stability robustness 
objectives are often difficult to achieve via the minimization of a quadratic performance index, and 
that the resulting control laws are often unnecessarily complex. LQG control laws are even less 
desirable in the multirate as compared to the single-rate case because multirate Kalman filter and 
LQR state feedback gains are periodically time-varying [2]. In short, LQG multirate sampled-data 
control laws can provide a useful benchmark for performance comparisions, but they are not 
practical for applications. 

Parameter optimization methods [2,19] for multirate sampled-data control law synthesis 
combine the principal advantages of the LQG and successive loop closures synthesis methods. They 
allow the synthesis of multirate sampled-data control laws of practical structure, and 
simultaneously account for all dynamic coupling between control loops. The typical parameter 
optimization method requires that the control law structure and its parameters to be optimized be 
prescribed. A numerical search is used to determine values for those parameters such that a 
performance index is minimized, possibly subject to constraints on those parameters. The 
disadvantage of parameter optimization methods is that they inevitably require a numerical search 
to determine the control law parameters. 

A new parameter optimization method for synthesizing multirate sampled-data control laws is 
described in Sec. Ill of this paper. It is the second generation of the method described in Refs. 2 and 
20. Unlike its predecessor, which accomodates only partial state feedback control laws, this new 
method accomodates a general, dynamic, multiple-input multiple-output control law structure. 
This new control lawstructure is described in Sec. II of this paper. Section IV describes an 
application of this new method to a design problem involving a two-link robot arm model. 
Conclusions are given in Sec. V. 
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II. CONTROL LAW STRUCTURE 


This section describes the multirate sampled-data control law structure in Fig. 2. In Fig. 2, y is 
the noise-free, continuous-time sensor signal, v is the discrete-time sensor noise signal, and u is 
the continuous- time control signal. The one sampler in Fig. 2 operates at the sampling rate 1/T, 
where T is the STP of the system's complete sampling/update schedule. The Delay blocks are one- 
STP delays. The ZOH block is a zero-order hold. 

The sensor samplc-and-hold dynamics are represented by 

y(m,n+l) = II - Sy(n)] y(m,n) + Sy(n) y(m,n) ^ 

where y is the sensor signal hold state vector. The matrix Sy(n) is the sensor switching matrix for 
the (n+ l)th STP. We define a switching matrix as a diagonal matrix with 1 or 0 at every diagonal 
position. If the ith diagonal element of Sy(n) is 1, the continuous-time signal from the ith sensor is 
sampled at the start of the (n+ l)th STP of every BTP and that sampled value is immediately stored 
as the ith element of y; otherwise, the same element of y is held at those instants. The key point is 
that y always contains the most recent sampled sensor data. 

The processor dynamics are represented by 

z(m,n+l) = [I - S 7 (n)]z(m,n) + S z (n) (A z (n) z(m,n) 

+ B z (n) {[I - Sy(n))y(m,n) + Sy(n) y(m,n)j] (2) 

u(m,n) = C z (n) z(m,n) 

+ D z (n) ([I - Sy(n)] y(m,n) + Sy(n) y(m,n)} ^ ) 

where z is the processor state vector, and 0 is the processor output vector. The matrix S z (n) is the 
processor state switching matrix. If the ith diagonal element of S z (n) is 1, the ith processor state is 
updated at the start of the (n+l)th STP of every BTP; otherwise, the same element of z is held at 
those instants. The matrices A z (n), B z (n), C z (n), and D z (n) are the processor state model matrices. 
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whose determination constitutes the control law synthesis problem. Note that a nonzero D z (n) 
results in direct feedthrough of sensor data to u(m,n). 

The control signal update-and-hold dyamics arc represented by 

u(m,n+l) = [f - S^n)] u(m,n) + S^n) u (m,n) (4) 

where u is the control signal hold state vector. The matrix S^n) is the control signal switching 
matrix. If the ith diagonal element of Su(n) is 1, the ith element of u is updated at the start of the 
(n+ l)th STP of every BTP; otherwise the same element of u is held at those instants. 

Finally, the continuous-time control signal u is generated by 

u(t) = [I - S^n)] u(m,n) + S u (n) u(m,n) (5) 

for all t on [(mP + n)T, (mP + n + 1)T). 

The advantage of the control law structure of (1) through (5) is that it can be used to represent 
virtually any sampled-data control law structure of practical interest Its form, however, is not 
standard. Straightforward algebra, applied to (1) through (5), yields the following more standard 
form: 

c(m,n+l) = A c (n) c{m,n) + B c (n) y(m,n) (6) 

u(m,n) = C c (n) c(m,n) + D c (n) y(m,n) (7) 

where 

c(m,n) = [z(m,n) u(m,n) y(m,n)l T (8) 
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(9) 


A c (n) = 


[I - S z (n)l + S z (n) A z (n) S z (n) B z (n) (I - S y (n)| 


0 


I - S y (n) 


Su(n) C z (n) S u (n) D z (n) [I - Sy(n)| 


0 

0 

I - S u tn) 


B c (n) = 


S z (n) B z (n) Sy(n) 
S y (n) 

|_S u (n) D z (n) Sy(n)J 


C c (n) = (S^n) C z (n) S u (n) D z (n) [I - Sy(n)] I - S^n)] 


D c (n) = [SJn) D z (n) S y (n)] 


with 


(10) 


( 11 ) 

( 12 ) 


u(t) = u(m,n) 

(13 

for all t on ((mP + n)T, (mP + n + 1)T). 

III. PARAMETER OPTIMIZATION METHOD 

This section describes a parameter optimization control law synthesis method for the control 
law structure of Sec. II. It is a generalization of the similar method for state feedback control laws 
described in Refs. 2 and 20, and incorporates also the multiple-plant-condition design for 
robustness ideas of Ref. 21. The approach involves a numerical search to determine the processor 
matrices, A z (n), B z (n), C z (n), and D z z(n), for n=0, . . . , P-1, such that a quadratic performance index 
is minimized. That approach has been criticized in the past because of (1) the difficulties of 
achieving practical performance and stability robustness objectives via the minimization of a 
quadratic performance index, and (2) difficulties related to the convergence of the numerical search. 


The proposed method addresses those criticisms in several ways. First, to enable synthesis for 
robustness to plant parameter variations, the performance index is defined over multiple plant 
conditions. This simple idea has been a key to the success of the popular Sandy [21,22,23,24,25,26] 
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algorithm for synthesizing robust continuous-time control laws. Second, to improve the 
convergence of the numerical search, the performance index and its gradients with respect to the 
control law parameters are calculated exactly, at every iteration, using closed-form expressions. 
Third, so that a stabilizing initial guess for the control law is not required, and to eliminate 
problems with destabilizing control laws encountered during the search, a finite-time performance 
index is used. Finally, to lessen the difficulties of achieving practical performance and stability 
robustness objectives via the minimization of a quadratic performance index, linear and nonlinear 
constraints can be imposed on the control law parameters. 

The continuous-time plant dynamics at plant condition i are assumed to be represented by. 


(t) = Ap’p'^l) + S“u ,il (l) + B^* 0 > (t) 


(14) 


y (i) (t) = p (i) (t) 


(15) 


where p (i) is the plant state vector, u (i) is the control input vector, y (,) is the sensor output vector, 
and w (i) is a stationary, zero mean, gaussian white noise input vector of known power spectral 
density. 


The performance index is assumed to be 


J = 



Q (i> 0 
0 R (i) 


p (i) (t) 

u (i) (t) 


dth 


(16) 


where N p is the number of plant conditions; E is the expected value operator; t f is the final time 
and is a multiple of the BTP of the system's complete sampling/update schedule; and Q (l> and R (l) 
are the state and control weighting matrices for the ith plant condition and are non-negative 

definite matrices. 
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Based upon the description of the continuous-time plant dynamics in (14) and (15), a complete 
description of the complete system s sampling/ update schedule, the performance index in (16), and 
the control law in (6) through (13), closed-form expressions for the performance index J and for its 
gradients with respect to the processor matrices A z (n), B z (n), C z (n), and D z (n), for n=0, . . . , P-1, arc 
derived in Refs. 19 and 27. Those derivations and the resulting closed-form expressions are 
lengthy, and will not be repeated here. The key points are that the resulting expressions are closed- 
form, and that the number of computations required for their evaluation is independent of t f . The 
single restriction for those expressions to be valid is that the state transition matrix for the BTP-to- 
BTP closed-loop system must be diagonalizable [19]. That is not a serious restriction because that 
matrix is rarely nondiagonalizable in practice. 

Thus far nothing has been said about synthesizing other than periodically time-varying control 
laws. To that end, the performance index and gradient derivations in Refs. 19 and 27 assume that 
the processor matrices are constrained to satisfy 


' D z (n) C z (n) 

M-l 

D z (r) C z (r) " 

. B z (n) A z (n) 

= £ a(n,r) 

r=0 

_ B z (r) A z (r) . 


with M e { 1 , . . . , p), and with the a functions constrained to satisfy 

a(n,p)a(n,q)={’ (18) 

hquations (17) and (18) constrain the number of different sets of processor matrices to M The 
function a(r,n) determines which set of processor matrices is active at the (n+ l)th STP. Equation 
(18) guarantees that only one set of processor matrices is active per STP. 


Based on the description of the continuous-time plant dynamics in (14) and (15), a complete 
description of the complete system's sampling/update schedule, the performance index in (16), the 
control law in (6) through (13), the constraint relations in (17) and (18), and the closed-form 
expressions for the performance index J, and for its gradients with respect to the processor matrices 
A z (r), B z (r), C z (r), and D z (r), for r=0, . . . , M— 1, in Refs. 19 and 27, we have developed a computer 
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algorithm to numerically determine a set of processor matrices that minimizes J. A numerical 

search is used to determine the processor matrices A z (r), B z (r), C z (r), and D z (r), for r = 0, M-l, 

given an initial guess for those matrices. The NPSOL nonlinear programming algorithm is used 
for the numerical search. NPSOL [28] is a powerful nonlinear programming package with good 
convergence properties as a result of its use of exact performance index and gradient evaluations at 
every iteration. In addition, NPSOL accomodates linear and/or nonlinear constraints on the 
independent variables. This means that linear and/or nonlinear constraints on the control law 
parameters can be combined with the usual performance index minimization objectives to achieve 
practical performance and stability robustness objectives. 

Additional important features of this new synthesis algorithm include automatic discretization 
of the continuous-time plant model and of the continuous-time performance index [27], These are 
important design features because they effectively decouple the sampling/update rates selection 
problem from the problem of determining a suitable performance index. This means that the 
performance index can be determined first, based on a continuous-time design, and that this new 
algorithm can then be used to determine a multirate sampled-data design that minimizes the same 
performance index. 

In summary, the inputs required to apply this new synthesis algorithm are the following: 

• A state model description of the continuous-time plant dynamics at each of the N,, 

plant conditions. 

• State and control weighting matrices for the performance index at each of the N p plant 

conditions. 

• The final time tf for the performance index. 

• The power spectral density of the continuous-time white process noise at each of the 

Np plant conditions. 

• A complete description of the complete system's sampling/ update schedule. 
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• The integer M and the a functions that constrain the periodicity of the processor 

matrices via (17) and (18). 

• The desired dynamic order and structure for the processor matrices. 

• The covariance matrix for the discrete-time sensor noise at each of the N p plant 

conditions. 

• A complete description of all linear and/or nonlinear constraints to be imposed on the 

elements of the processor matrices. 

• An initial guess for the processor matrices. 

A disadvantage of most parameter optimization control law synthesis methods is that they 
require a stabilizing initial guess for the control law. That is not the case with this method because 
of its finite-time performance index. The finite time ensures that the performance index and its 
gradients will be finite whether or not the closed-loop system is stable. A disadvantage of the finite 
time is that a steady-state solution, i.e., for 

Jss= * im J (19) 

cannot be obtained directly. A steady-state solution is easily obtained, in practice, however, by 
choosing a finite time t f that is large compared to the characteristic times of all of the closed-loop 
system’s poles. Because the number of computations required to evaluate the performance index 
and gradient expressions of Refs. 19 and 27 does not depend upon t f , this can be done without 
penalty in terms of the computation time for the numerical search. 

In practice, because digital computers cannot store arbitrarily large finite numbers, a steady-state 
solution usually cannot be obtained by simply initially setting tf to a large value. Instead, it is 
usually necessary to complete first (i.e., when the current best guess for the control law parameters 
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is poor) an optimization for a small tf, and to then re-optimize, for larger and larger tf, until t f gets 
to be large compared to the characteristic times of all closed-loop poles. 


The final key issue regarding this new synthesis algorithm concerns the requirement that the 
structure of the processor matrices be specified. The key point is that the imposed structure should 
guarantee that the free parameters for the numerical search constitute an independent set with 
respect to the control law's input-output dynamics. When the processor dynamics in (2) and (3) are 
considered, with the constraints in (17) and (18) in effect, it is straightforward to see that the 
complete set of the elements of A z (r), B z (r), C z (r) and D z (r), for r = 0, . . . , M-l, do not constitute 
such an independent set because, for example, an arbitrary change in one element of B z (0) can be 
compensated for by changes to the elemets of C z (0), and to the other elements of B z (0), such that the 
processor's input-output dynamics are unchanged. 

Thus, additional structure, or, equivalently, additional constraints, must be imposed on the 
elements of the processor matrices to guarantee that the free parameters for the numerical search 
constitute an independent set with respect to the control law's input-output dynamics. In practice, 
a suitable set of such constraints can frequently be determined based on "classical" control law 
structures (e.g., combinations of lead and lag compensators and notch filters). 


More generally structured control law can, of course, also be accomodated. What constitutes an 
optimal structure for the processor matrices for the general case is a topic of current research. We 
have successfully applied the following structure (shown for the n-is-even case) for the particular 
case where the constraints in (17) and (18) are applied with M = 1 (the time-invariant case): 

T 0 1 “1 


A z (0) = block diag 




-20i 


f, i = 1, . . . , n/2 


( 20 ) 


B z (0) = 


b n 

bnl 


b|m 

bnm _ 


( 21 ) 
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C 7 (0) = 


1 

C 21 
L Cp] 


C2n 


l -pm -i 


DJO) = 


*!1 


L d pl 


d 


pm J 


( 22 ) 


(23) 


The appendix shows that the (m+ p)n+ mp a„ 0 i { , bjj, <- and d sj parameters of this structure 
constitute an independent set with respect to the control law's input-output dynamics (for any 
non-trivial sampling/update schedule) provided that no e.genvalues a* ± jeoj of the processor 
dynamics are repeated. 


IV. EXAMPLES 


This section describes the design of a tip position control system for a planar two-link robot 
arm. The robot arm system is shown in Fig. 3. The first link is long and massive, for large-scale 
slewing motions. The second is relatively short and lightweight, so that high-bandwidth control of 
the arm's tip position can be achieved using a relatively small motor at the second joint. The pin 
joint, rotational spring and rotational damper at the midpoint of the first link models flexibility in 
that link. The second link is assumed to be rigid. The motor torques T, and T 2 are the control 
inputs, and it is assumed that only the joint angle 0 and the tip position 5 are measured. The 
linearized dynamical equations for this system for small e - 0 and small <j> - e are easily derived. 

The spring constant (k) and damping coefficient (b) values (in Fig. 3) were chosen based on that 

model to achieve 1 percent damping and a 10 Hz natural frequency for the open-loop vibration 
mode. 

Figures 4 through 6 show the closed-loop arm responses, based on the linearized arm 
dynamics, to a step change in the commanded tip position with nine different control laws. The tip 
position (5) responses are shown in Fig. 4. The simultaneous control torque (T, and T 2 ) responses 
arc shown in Figs. 5 and 6. The nine different control laws are briefly described as follows: 
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LQR Analog: The continuous-time LQR (full-state-feedback) control law that minimizes 


}= lim 2f 

tf ~ > « ‘ 


2[a2 , 


(^[6 + (p5) 1 + 


^ 1 Imax j 


t 2 v 


[ 2 max , 


dt 


( 24 ) 


with 


p = 32 rad -1 


(25) 


p = 14 rad/m 


(26) 


T 2 max = 0.00335 N m 


(27) 


Tlmax/T^max - 8 


(28) 


The T 2max value is the T 2 torque that achieves^ = 2 n rad/sec in 1 sec with 0(t) s e(t) s 0 and<}i(0) = 0. 
The T lmax /T 2max value represents a typical ratio of peak motor torques at the respective joints. The 
p and p values were chosen by trial and error to achieve the closed-loop poles in Table 1. Note that 
the ratio of the characteristic frequencies of the rigid-body closed-loop pole pairs is eight; and that 
the characteristic frequency of the faster closed-loop rigid-body pole pair is a factor of five less than 
the characteristic frequency of the closed-loop vibration mode. 

Third-Order Analog Successive Loop Closures: The third-order, continuous-time, successive loop 
closures control law in Fig. 7, which consists of a single lead compensator in the 0-to-T, loop and 
twin, cascaded lead compensators in the S-to-T 2 loop. The closed-loop poles for this design are in 
Table 2. Note that the rigid body and vibration mode closed-loop poles match those of the LQR 
Analog design. 


Third-Order Multirate Tustin: A multirate sampled-data approximation to the Third-Order 
Analog Successive Loop Closures design obtained via Tustin’s approximations of the continuous- 
time transfer functions in Fig. 7. The sampling/update rates (in samples/updates per second) of the 
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e-to-T, and 5-to-T 2 loops are eight times the characteristic frequencies (in cycles per second) of the 
slow and fast, respectively, rigid body closed-loop pole pairs from the Third-Order Analog 
Successive Loop Closures design. 


Optimized Third-Order Multirate Tustin: The same as the Third-Order Multirate Tustin design, 
but with the lead compensator gain, zero, and pole locations optimized, by the paramether 
optimization control law synthesis method of Sec. Ill, to minimize the same performance index as 
in the LQR Analog design. To synthesize this control law, continuous-time process noise and 
discrete-time sensor noise inputs were added to the robot arm model. The former were taken to be 
white noise disturbance torques w, and w 2 , coincident with the respective control torques, with 



(w^t) w 2 (t)| ) = 


4.9x1 0" 5 0 

0 1.6xl0 -5 


5(t-x) 


(29) 


where 5 is the Dirac delta function. The latter were taken to be stationary, purely random 
sequences, V} and v 2 , for the 0 and 6 measurements, respectively, with 



v^mm) " 
v 2 (m,n) . 


[Vj(m,n) v 2 (m,n)J } = 


8.1 xlO -5 0 
0 lxlfT 4 


(30) 


Multirate Third-Order: The same as the Optimized Third-Order Multirate Tustin design, but using 
the third-order, generalized, time-invariant structure in (20) through (23) for the processor 
matrices. Just as with the Optimized Third-Order Multirate Tustin design, two of the processor 
states are updated at the faster sampling/update rate, and the third is updated at the slower 
sampling/ update rate. 


Multirate Second-Order: The same as the Multirate Third-Order design, but using the second- 
order, generalized, time-invariant structure in (20) through (23) for the processor matrices. One of 
the processor states is updated at the faster sampling/update rate, and the other is updated at the 
slower sampling/update rate. 


43 



Multirate First-Order: The same as the Multirate Third-Order design, but using the first-order, 
generalized, time-invariant structure of (20) through (23) for the processor matrices. The one 
processor state is updated at the faster sampling/ update rate. 

Single-Rate Third Order: The same as the Multirate Third-Order design, but single-rate, with the 
single sampling/update rate chosen to yield the same number of real-time computations per unit 
time as the Multirate Third-Order design. 

Analog Third-Order: The continuous-time equivalent to the Multirate and Single-Rate Third- 
Order designs. The processor matrices have the same structure as in the Multirate and Single-Rate 
Third-Order designs. The control law was synthesized using the Sandy algorithm [21] to minimize 
the same performance index as in the Multirate and Single-Rate Third-Order designs. 

The LQR Analog responses in Figs. 4a, 5a and 6a constitute the optimal responses for the 
performance index in (24), assuming full state feedback, no process or sensor noise, and infinitely 
fast sampling. The Third-Order Analog Successive Loop Closures responses in the same figures 
have low tip position overshoot, but include also a relatively large contribution from the vibration 
mode (see especially Figs. 5a and 6a). 

The Third-Order Multirate Tustin responses in Figs. 4a, 5a and 6 are unacceptable. This is 
somewhat suprising, but not totally unexpected given the (low) factor-of-eight sampling/update 
rate-to-characteristic frequency ratio for this design. 

The Optimized Third-Order Multirate Tustin responses in Figs. 4a, 5a and 6a are acceptable, and 
demonstrate that the parameter optimization control law synthesis algorithm of Sec. Ill can be used 
to optimize the parameters of classically-structured control laws. 

The Multirate Third-Order, Second-Order and First-Order responses in Figs. 4b, 5b and 6b 
demonstrate that the same parameter optimization control law synthesis algorithm can be used to 
synthesize multirate sampled-data control laws having a prescribed dynamic order and a 
prescribed, but general, structure, with apriori specified sampling/update rates for all sensors, 
processor states, and control inputs. 
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The Single-Rate Third-Order and Analog Third-Order responses in the same figures put the 
multirate responses in perspective. The Single-Rate Third-Order control law is the single-rate 
equivalent to the Multirate Third-Order control law because it (1) was synthesized to minimize the 
same performance index, using the same process and sensor noise characteristics; and (2) requires 
the same number of computations per unit time for real-time operation. 

The Analog Third-Order responses in the same figures are the responses that would have been 
obtained with the Multirate Third-Order control law and Singe-Rate Third-Order control laws if 

sampling and update rates were not an issue, and very fast sampling and update rates were 
everywhere used. 

V. CONCLUSIONS 


With the possible exception of successive loop closures, the multirate sampled-data control law 
synthesis methods available today fail to provide the designer with sufficient flexibility to prescribe 
sensor sampling rates and processor state and control input update rates. A new parameter- 
optimization-based method for synthesizing multirate sampled-data control laws of arbitrary 
dynamic order that provides that flexibility is described in this paper. This new method, described 
in Sec. Ill, determines, by numerical optimization, the free parameters of the general purpose 
multiple-input, multiple-output, sampled-data control law structure in Fig. 2, to minimize a 
quadratic performance index, possibly subject to linear and/or nonlinear constraints on those 
parameters. A stabilizing initial guess for the control law is not required because the performance 
index is finite-time. To enable the synthesis of robust control laws, the performance index can be 
defined over multiple plant conditions. 

An application of this new method to the design of a tip position control system for a sixth- 
order, two-link robot arm was described. Multirate sampled-data control laws of various dynamic 
orders synthesized by various methods were compared to confirm that the new synthesis method 
can be used to synthesize multirate sampled-data control laws having a prescribed dynamic order 

and structure, with apriori specified sampling/update rates for all sensors, processor states, and 
control inputs. 
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APPENDIX A. SUPPORTING MATHEMATICS 


Consider the control law in (1) through (5). Suppose that the constraints in (17) and (18), with 
M=l, are in effect, so that the processor matrices are constrained to be time-invariant. Suppose that 
the processor matrices, A z (0), B z (0), C z (0) and D z (0), are further constrained to have the forms in 
(20) through (23). Finally, to guarantee a nontrivial sampling/update schedule, suppose that the 
sensor, processor state, and actuator switching matrices satisfy 


del 


P-1 

IV n) 

n=0 


* 0 


det 


P-1 

£s z (n) 

n=0 


* 0 


det 


P-l 

£s u (n) 

n=0 


* 0 


(31) 


(32) 


(33) 


We will show that the (m+p)n+pm a,, cd , b;j, cjj and djj elements of the control law then constitute 
an independent set with respect to that control law's input-output dynamics if and only if A z (0) has 
no repeated eigenvalues. 

We begin by noting that, with (31), (32) and (33) in effect, it is straightforward to see that the 
independence in question does not depend whatsoever on the sensor, processor state, or actuator 
switching matrices. Therefore we consider only the special case where S y (n), S z (n) and S u (n), for 
n=0, . . . P-l, are identity matrices. The control law then reduces to 

z(m,n+l) = A z (0) z(m,n) + B z (0) y(m,n) (34) 

u(m,n) = C z (0) z(m,n) + D z (0) y(m,n) (35) 


where u(m,n) is defined in (13). 


48 



Consider first the control law 


z(k+l) = A z(k) + By(k) 


136 ) 


u(k) = C z(k) + D y(k) 


(37) 


where 


z = 


"zf 


~yf 


" u f 


- z n- 

y = 

-ym- 

U - 

- u p- 

(38) 


A = 


0 

0 


0 


(39) 


and B, C and D have the forms in (21) through (23). So that the control law’s impulse response will 
be purely real, the X-s and the corresponding columns of C and rows of B must be either real, or 
must occur in complex-conjugate pairs, and D must be purely real. 

Lemma: The (m+p)n+pm X u b iy Cjj, and djj parameters (counting a real element as one parameter, 

and a complex conjugate pair of elements as two parameters) of the control law in (36) and (37) 

constitute an independent set with respect to that control law’s input-output dynamics if and only 
if for i*j. 

Proof: Consider the related control law 


z(k+l) = A z(k) + B y(k) 


(40) 


u(k) = Cz(k) 


( 41 ) 


with 
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(42) 


C = 


<=u 


c P l 


c In 
'-pn - 


Its input-output dynamics arc represented by 


n 


H(z) = C (z I-Ar 1 



CjBj 

z-Xj 


(43) 


where C t is the i th column of C, and B* is the ith row of B. The bij and Cjj parameters of this 
control law are dependent with respect the control law's input-output dynamics if and only if, for 
an arbitrary change in one, the others can be changed so that H(z) is unchanged. 


Case 1: No repeated Xj ’s. 

From (41), for the case of no repeated X, 's, it is straightforward to see that when one of the X; s is 
changed by an arbitrary amount, it will not be possible to change the remaining Xj, b,j and c,j 
elements so that H(z) is unchanged. But when bjj is multiplied by a nonzero but otherwise 
arbitrary a, we can multiply the remaining elements of B[ by a, and divide Cj by a, so that H(z) is 
unchanged. Thus, the Xj, bjj and c S j parameters of the control law in (40) and (41) are dependent 
with respect to that control law's input-output dynamics. 

If, however, one element of every column of C is fixed, as is the case in the C matrix of (22), it is 
similarly straightforward to see that it will not be possible to compensate for an arbitrary change in 
any Xj, bjj or c t j element by changing the remaining Xj, bjj and Cjj elements so that H(z) is 
unchanged. Thus, for the case of no repeated Xj’s, with one element of every column of C fixed, 
the Xj, bjj and Cjj parameters of the control law in (40) and (41) constitute an independent set with 
respect to that control law’s input-output dynamics. 


Case 2: Repeated Xj's. 

Consider again the control law in (40) and (41), but this time suppose that Xj=X 2 . That control law's 
input-output dynamics are represented by 
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(44) 


H(z) = C (z I-A)- 1 


B = 


Cj B 1+ C 2 B2 
Z-Xj 


n 



Cj Bj 
z-Xj 


From (42), it is straightforward to see that, with or without one element of every column of C fixed, 
the remaining bjj and Cjj elements of Cj, Bj,C 2 , and B 2 can be changed to compensate for an 
arbitrary change in any one element of C v Bj, C 2 , or B 2 so that H(z) is unchanged. Thus, for the 
case of repeated Xj's, with or without one element of every column of C fixed, the Xj, b^ and cjj 

parameters of the control law of (40) and (41) are dependent with respect to that control law s input- 
output dynamics. 


General Case: We conclude that the (m+p)n Xj, bjj, and qj parameters of the control law in (36) and 
(37), with D= 0, constitute an independent set with respect to that control law s input-output 
dynamics if and only if Xj^Xj, for i/j. A nonzero D matrix simply adds pm parameters to that set. 

Theorem: For the control law in (1) through (5); with the constraints in (17) and (18), with M=l, in 
effect, so that the processor matrices are constrained to be time-invariant; with A z (0), B z (0), C z (0) 
and D z (0) further constrained to have the forms in (20) through (23); and assuming that the sensor, 
processor state, and actuator switching matrices satisfy (31) through (33); the (m+p)n+pm a ; , cOi, b Sj , 
c^ and djj elements of that control law constitute an independent set with respect to that control 
law s input-output dynamics if and only if A z (0) has no repeated eigenvalues. 


Proof. The control law in (34) and (35) has the same number of free parameters as the control law 
in (36) and (37), and the two are related by the similarity transformation x= M z , where 


M = block diag { — 


CTi+jUj 1 
-Oj+jCOj -1 


i = 1, . . . , n/2 ) 


(45) 


with <Tj = RefXj) and C0j = ImKXj). 
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Closed- Loop 
Pole 

Damping 

Ratio 

Charactenstic 

Frequency 

Rigid Body 

-1.10 ±./‘ 1.10 

0.71 

0.25 Hz 

Rigid Body 

-8.81 ±j 8.83 

0.71 

2.0 Hz 

Vibration Mode 

-0.649 ±j 62. 8 

0.01 

10 Hz 


Table 1 LQR Analog Design Closed-Loop Poles 



Closed-Loop 

Pole 

Damping 

Ratio 

Characteristic 

Frequency 

Rigid Body 

-l.10iyl.il 

0.71 

0.25 Hz 

Rigid Body 

-8.88 ±j 8.84 

0.71 

2.0 Hz 

Vibration Mode 

-1.35 ±j 63.9 

0.02 

10 Hz 

Compensator 

-10.5 

- 

1.7 Hz 

Compensator 

-33.2 ±j 34.0 

0.70 

7.6 Hz 


Table 2 Third-Order Analog Successive Loop Closures Design Closed-Loop Poles 
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Time Lines for Sampling/Update Activities: 



4T 8T 12T 16T 20T 24T 



0 3T 6T 9T 12T 15T 18T 2 IT 24T 


BTP 


Fig. 1 Example Multirate Sampling/Update Schedule 
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Parameters: Mass Length 

L, 0.5 kg 0.5 m 

L 2 0.5 kg 0.5 m k = 37.33 N/rad 

L 3 0.04 kg 0.2 m b = 0.012 N- s/m 

The natural frequency of the vibration mode is 10 hz. 

Inputs: Torques T x and T 2 

Outputs: 0 and 8 


Fig. 3 Two- Link Robot Arm 
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LQR Analog 

Third-Order Analog Successive Loop Closures 
Third-Order Muliirate Tustin 
Optimized Third-Order Multirate Tustin 


Seconds 


Fig. 4a Tip Position (5) Responses to Tip Position Step Command 


-Multirete Second-Order 
— Multintfa Third-Order 


Analog Third-Order 

Multirate First-Order 

Multirate Second-Order 

Multirate Third-Order 

Single-Rate Third-Order 


Fig. 4b Tip Position (8) Responses to Tip Position Step Command 






Netwon-Meters Netwon-Meters 
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Netwon-Meters 
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Fig. 7 Third-Order Analog Successive Loop Closures Design 
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Reduced Order Multirate Compensator Synthesis 

Gregory S, Mason* and Martin C. Berg** 

University of Washington 


Abstract 

A method for synthesizing reduced order mullirate compensators is presented. Necessary conditions for 
which the compensator parameter values minimize an infinite lime quadratic cost function are derived. An 
algorithm for finding compensator parameter values which satisfy the necessary conditions is described This 
algorithm is then used to design several tip position controllers for a two link robot arm. 

INTRODUCTION 

In many cases a multirate compensator can provide better performance than a single rate compensator 
requiring the same number of real-time computations. Berg, for example, was able to reduce the steady state 
RMS response of states and controls to a disturbance for a simple mass-spring-mass system nearly 20% by 
using a multirate compensator over a single rate compensator [1]. Numerous other examples have been 
provided in the literature by Berg [1 ]-[3], Amit [4]-[5], and Yang [6]. While multirate compensators can 
provide improved performance over single rate compensators, they are also, in general, more complicated to 
design. 

The complexity of multirate compensators stems from the fact that they are by nature time varying, 
periodically time varying for most practical applications. Not only must designers choose multiple 
sampling/update rates for the compensator, but they must also determine the parameter values for a time 
varying compensator. 

One method for designing multirate compensators is multirate LQG [4], Multirate LQG is the multirate 
equivalent of single rate LQG and is straightforward to solve because the equations governing the solution are 
similar to those for the single rate case. Multirate LQG, however, results in a full order compensator which 


* Graduate Student, Department of Mechanical Engineering FU-10, University of Washington. Seattle WA 98195 

* * A1AA Member, Assistant Professor, Department of Mechanical Engineering FU-10, University of Washington, Seattle 

WA 98195 
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has periodically time varying gains. For many applications full order time varying compensators arc not 
practical. 

A Generalized Algorithm for Multirate Synthesis (GAMS) [6| was developed by Yang to overcome many 
of the shortcomings of multirate LQG. Yang’s algorithm can synthesize reduced order multirate compensators 
with or without time varying gains by using a numerical gradient type search to find optimum compensator 
parameter values. His algorithm uses a finite time cost function in its problem formulation, unlike multi rate 
and single rale LQG which use an infinite time cost function. By using a finite time cost function Yang’s 
algorithm eliminates the numerical problem that arises when a destabilizing compensator is encountered during 
the numerical search. Even though Yang’s algorithm uses a closed form expression for the gradient, the 
calculations necessary to perform the gradient-type search are extremely cumbersome. 

In this paper we present a new algorithm for synthesizing reduced order multirate compensators with or 
without lime varying gains. The algorithm utilizes the compensator structure of Yang’s algorithm, but the 
problem is formulated using an infinite lime, instead of a finite time, cost function. This allows us to derive 
necessary conditions for which the multirate compensator minimizes the cost function. The equations for the 
necessary conditions are fairly simple and can be solved directly using a standard nonlinear equation solver, 
eliminating many of the numerical complexities of Yang’s algorithm. 

The General Multirate Compensator 

Before deriving the equations governing a reduced order multirate compensator, we will first present the 
structure for a general multirate compensator. We restrict our discussion for now to compensators with time 
invariant gains and sampling/update rates whose ratios are rational numbers. 

A general multirate compensator is shown in Figure 1. Each input (y), output ( u ), and state (z) is 
sampled/updated at a rate which, in general, represents the desired bandwidth of the input or output with 
which it is associated, y is the value of y currently available to the digital processor from the zero order hold; 
while u is the current output from the digital processor which is held with a zero order hold to form the output 
u. When the sampling/update rates have ratios which are rational numbers the sampling/update schedule is 
periodically time varying. We define the greatest common divisor of all the sampling/update periods as the 
.shortest time period (STP) and the least common multiple of all the sampling/update periods as the basic time 
period (BTP) (see Figure 2). 
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The state equations for the multirate compensator pictured in Figure 1 are: 


) 

[I-Sz.k]+S z ,kA 

S 2tk B[l-Sy X ] 

0 

iz 1 

$z,kB Sy.k 

- 

0 

[I-Syjt] 

0 

\l\ + 

Sy,k 


$u,kC 

Su t kD[I~Sy,k\ 


\ U Ik 

$u,k&S y,k 


u k - [s u ,kC s Ut kD[!’Sy,k] [v*^*,*] y* (2) 


u is a hold slate used to model the sampler and zero order hold between u and “• V’ and v* are 
switching matrices for y f z> and u respectively that model the system’s sampling/update activity at the start of 
the k th STP. s+ k has the form: 


s *,k = 


r x 0 0 ... 0 

0 r 2 0 ... 0 

0 ... 0 r m ^ x 0 

0 . 0 0 r m * 


f 1 if the j th (z\ y, or u) is sampled/updated 
where r j = \ at the start of the k tfl STP 

lO otherwise 


m z = the number of states (z) 
rriy - the number of inputs (y) 
m u - the number of outputs (u) 


A more complete discussion of this compensator structure can be found in [6]-[7]. 
Equations (1) and (2) can be written more compactly as: 


z k + \ = A k z k + B k>k o> 

u k = c \?k + D k>k (4) 
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where z k = 



Equations (3) and (4) form a single rate periodically time varying system with a sampling rate of one STP 
and a period of one BTP. If jV=BTP/STP, then A k = A^ +iV , B k = B k + N , C k = and D k = D k ^. 

Even though A k , B k , C k and D k are periodically time varying, the multirate compensator gains, A, E , C> 
and D , are time invariant. The periodicity of the multirate compensator is due to multirate sampiing/updating, 
not the compensator gains. In the remainder of this section we will demonstrate how the time invariant 
compensator gains, A t 5, C, and D, can be separated from the periodic compensator matrices A k , B k , C k and 

% 

Define the composite compensator matrix as 


Pk = 


D k 

B k 


C k 

A k 


(5) 


and factor P k as follows: 


where 


P k ~ 5 1^2* + S 3k 


P = 


D C' 

-B A. 


S\ k 


Sujl 0 

0 Sz'k 

0 0 

$u,k 0 


■S' 2* = 


Sy,k 

0 


0 I-Syjc 0 
/ 0 0 


S3* = 


0 0 0 l-s uJc 

0 I -s 2tk 0 0 

Sy,k 0 /-Syj[ 0 

0 0 0 l-Sujt 


(6) 

(7) 


( 8 ) 


(9) 


( 10 ) 


Equation (6) is a key result. It allows us to factor the time invariant compensator gains, the unknown 
parameters we will solve for in the next section, out of the time varying compensator. 

It is important to note the difference between P k and P in (6). P k (with a subscript) is a periodically time 

varying matrix defined by (5). It includes all the information about the compensator gains and the 
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.sampling/update schedule. P is a constant matrix which contains only the gains for the compensator. P k can 
be written in terms of P and S lk , S 2k , and S Jk using equation (6). S u , S 2k , and S^ k are periodically time 
varying matrices which contain a description of the sampling/update scheme. 

Derivation of the necessary Conditions 

In this section we will use the results of the previous section to derive the necessary conditions for the 
reduced order multirate compensator. The multirate problem to be solved is as follows: 

Given: the discretized plant model 


Xk+ 1 = Fxk + G«* + Ww k 


( 11 ) 


y k = Hx k + v* 


( 12 ) 


where F, G, W and H are obtained by discretizing the analog plant matrices at one STP; w k and v k are 

discrete-time Gaussian white noise inputs; u is the control input from the compensator, and y is the sampled 
sensor output. 

Find: the multirate control law with a prescribed dynamic order and sampling schedule, of the form of ( 1 )- 
(2), which minimizes a quadratic cost function of the form: 


'-s4'{er 


Q\ M 

M t q 2 


(13) 


where E is the expected value operator, and the summation from 1 to N accounts for the fact that the closed 
loop system is periodically time varying. A prescribed sampling schedule implies that the values of s 2 k , s k , 
and s u k are known. 

Using (3)-(4) it is easy to see that this problem is essentially a time varying feedback problem - a time 
invariant plant with a periodically lime varying compensator. One thing that makes this problem difficult is 
that the compensator has an explicit form, that of (l)-(2), in which only certain parameters, A, S, C, and D, 
can be adjusted to minimize J. 
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To solve the multirate control problem we cast it into output feedback form and follow a derivation similar 
to Mukhopadhyay’s for the single rate case [8]-[9]. Using (3)-(4), and (11 )-( 1 2) we write the output feedback 
equations: 


i i 


F 0 


G 0 


\ + 

( 

-0 0- 

1 z k | 

L o / Jl 



w* 

v* 


jy k 

\ Zk 


✓V 

H 

0 


0 

/ 1 

1 

- 0 

/ - 

W 

. 0 

oJ 

1 v* j 


I “* Lr o k c k ipd 

\ Z* +) j lB k A * J\ Z* / 


Equations (14)-(16) can be written more compactly as 


X M = Fx k + Gu k + Wr lk 

y k = Hx k + Vr ik 

u k = *V* 


(14) 


(15) 

(16) 


(17) 

(18) 
(19) 


It is important to keep in mind that P k in equation (19) corresponds to the P k in equation (5), a periodically 
time varying matrix which contains all the information about the multirate compensator gains and 
sampling/update rates. 


The closed loop system is 


x k+ i = F ck x k + G ck r] k 

(20) 

where 


F c k = F + GP k H 

(21) 

G ck =W+GP k H 

(22) 

The state covariance propagation for this system obeys: 


**♦! = F ck X k F T ck + G ck RG T ck 

(23) 
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where 


X k = E[x k x T k \ 

* = £{W) 


Equations (20)-(22) represent a periodically time varying system with a period of one BTP. We can 
generate a single rate system by repeated application of equation (20) over one BTP [10]. The single rate 
system can be written as 

*k+N = Fbk*k + Gbklbk (24) 


where 

Fbk = F c{k*N - 1 }Fc(k+N- 2)Fc(k+N-3y • -F c k 
Gbk = [Fc(k+N - 1 )F c(k+N-2Y ' -F c {k+\)Gck \ 

Fc(k+N-l)Fc{k+N-2yFc(k+2)Gc<,k+\)\ | 


(25) 

(26) 


T]bk ~~ 


rik + 1 




This single rate system has exactly the same values for* as the periodically time varying closed loop 
system at each BTP. However, the values of x at the intermediate STP’s are lost because x is incremented by 
N in (24) but only by 1 in (20). There are N such single rate systems associated with (20). They can be 
written as: 


Xk+u+i = F b{k +i>x k +i + Gb( k +i)r) k +i for i=l,2 ... >N 


(27) 


If F ££ is stable, then the periodically time varying system (20) is stable [11], We can calculate the steady 
state covariance for x using the following Lyapunov equations: 


x* = F bk X k F T bk + GbkRbGlk for £=1,2 N 


r /? 0 ... 01 


R b = 


0 R 0 


0 0 0 /? 


(28) 
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Note that X k is periodic, that is it varies within one BTP, but from BTP to BTP X k = X k+N . Once we 


have calculated X k at any k using (28), we can use (23) to propagate it over the BTP. This eliminates the need 
to solve equation (28) N times. 

Now, using (23), (13), and the properties of the Trace (TV) operator we can write the cost function for the 
stabilized system as (sec (8|-[9|) 


/ = X Htei + Mp S { + (MPkHf + (P k H?Q 2 P k H)x k + (P k V) r Q 2 P k VR) 


Adjoin the covariance constraints (23) to the cost J using Lagrange multipliers. A, t, to obtain: 

N 

J = X Tr[[Qi + MP k H + (MP k H) T + {P k H) T Q 2 P k H]x k + (P k VfQ 2 P k VR 

*= i (30) 


+ Al+\[F ck X k F T ck + G ck RGj k - X*+i]) 


with Xj = Xjy +1 . 


Necessary conditions for minimum J are 


57 = 0, - = 0, and ~ = 0 

5X* aA* +l 3/> 

a 2 / 

In addition, must be positive definite for a minimum J . 

a ? 2 

Substituting (30) into (31) and replacing F* with ^ rom (6) we obtain: 


= 0 = fli + M/y/ + [MPtf f + (PkHfQiPkH + F c r *A* +1 F c * - A* 

dX* 

for £=1,2,...,// with A^ = A A+/ y 


3A* +1 


= 0 = F ck X k Fc k + G ck RGl k - X k+i 


for k = 1.2....JV with = X k 


^L=0=2^jl k {[Q 2 +G T A k+i G]p k [HX k H T +VRV T ] 
dP *=1 

+[M 7 '+G r A* +1 F]jV*// r }5 2 7 * 
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Equations (32)-(34) arc a set of coupled matrix equations. They make up necessary conditions for P, 
which is comprised of the multirale compensator gain matrices A , B, C. and D , in equation (l)-(2), to 
minimize the cost function 7. Values of A , B, C. and D, found by solving (32)-(34) can be substituted into 
(l)-(2), along with the definition of the sampling schedule, s zk , s uk , and s y k , to form the complete time 
varying multirate compensator. 

To ensure that the compensator gains satisfying (32)-(34) minimize 7, we should also check that the 
Hessian of 7 with respect to P is positive definite. Our present algorithm docs not calculate the Hessian 

explicitly, but uses an approximate value calculated by the numerical search algorithm discussed in the next 
section. 

Equations (32)-(34) were derived assuming time invariant compensator gains. We can easily derive the 
corresponding equations for periodically time varying gain s Let 

A=A k , B = B k , C = C k , and D = D k ( 35 ) 


with the restriction that A k +„ = A k , B k+N = B k , C k+N = C k , and D k + N = D k 
Define the composite periodically time varying compensator matrix: 


Pk = 


D k 

B k 


C k 

A k 


(36) 


Then replace P with P k in (30) and differentiate with respect to P k to obtain 
Jr =0= 5 1 7 * l[Q2+G T A k+l G]p k [HX k H T +VRV T ] 

°rjt 

+[A 1 T +G T A k+i F]x k H T )sJ k for k = 1,2... TV 


(37) 


Thus for every new set of compensator gains we obtain one new equation of the form of (37). 

Equations (32)-(34) are very similar to the single rate equations. In fact, if we set S 2k and S 2k so 

they correspond to a single rate system, and N= 1, we obtain the exact results derived by Mukhopadhyay for 
ihc single rale case (8). 
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IMPLEMENTATION 


In order to find a reduced order multirate compensator that minimizes the cost function 7, we need to solve 
(32)-(34) for the compeasator gaias P. A flow chart of the algorithm used to determine the compensator gains 
is shown in Figure 3. Using the prescribed sampling schedule the algorithm first discretizes the analog plant 
model, analog cost function, and analog process noise model. See [2] for a discussion of the relevant 
discretization procedures. Equations (32)-(34) are then solved for the compensator gains using a gradient type 
search in MATLAB [12]. We chose a gradient type search to solve (32)-(34) because it allows us to easily add 
constraints on the parameters values - simple equality constraints were used to find the optimized 
compensators in the next section. The equations necessary to solve for the Lagrange multipliers are located in 
the Appendix, (A.3)-(A.4). To ensure that the solution represents a minimum 7, the algorithm checks that the 
Hessian of 7 with respect to the free parameters in P is positive definite at the solution point. 

Because (32)-(34) are not valid when the closed loop system is unstable, the algorithm 1) must be 
provided with an initial stabilizing compensator, and 2) must result in a stabilizing compensator at every 
iteration. From our experience, finding an initial stabilizing compensator is generally not a problem. Many 
systems suitable for multirate control can be stabilized using successive loop closure with minimal cross 
coupling between the control loops. A stabilizing multirate compensator can then be obtained by discretizing 
the individual continuous control loops at the desired sampling rates. When there are no constraints on its 
structure, a stabilizing compensator can also be obtained using the boot strapping method of Boussard [13]- 
[14], For difficult multirate control problems, where a stabilizing compensator cannot be found using either of 
the preceeding two methods, one can always use Yang’s algorithm to find a stabilizing compensator and then 
to switch to our algorithm to complete the optimization. In our experience Yang’s algorithm usually converges 
to a stabilizing solution quickly - it is the optimization of the compensator parameters that is time consuming. 

To avoid the problem of destabilizing compensators during the iteration process we included a check in the 
algorithm which systematically reduces the step size to ensure that the compensator is stabilizing. Because the 
gradient of the cost function with respect to the compensator parameters becomes very large near the stability 
boundary, the algorithm is always forced away from a destabilizing solution as long as it never steps over the 
stability boundary into an unstable region. 

Even though our algorithm was programmed as an interpreted Matlab M-File we found that it still 
performed better than Yang’s algorithm which runs as compiled FORTRAN. The primary difference between 
the two algorithms is in the complexity of the expression for the gradient of 7 with respect to the compensator 


70 



parameters. Calculation of the gradient expression for Yang’s problem involves diagonalization of the closed 
loop system and evaluation of several matrix equations with nested summations. Compare equations (32M34) 
with equations (112-115) in [31 to see the difference in the complexity of the two gradient expressions. 


Two Link Robot Arm Example 


We used a math model of a planar two link robot arm (TLA) to demonstrate the capabilities of our 
algorithm. This is the same model used by Yang [6], and so we were able to verify our results by direct 
comparison. A diagram of the TLA is shown in Figure 4. 

The goal of our design was to control the tip position (5) of the arm via a multirate compensator. We used 
the following analog cost function and process noise covariance matrices from [6], 


1 

" 0.21 

0 

0 

0 1 



J = limE{x T 

0 

0 

0 

0 

x + u T 

0.01 

° 1« 

l-*oo | 

0 

0 

18.5 

0 


0 

0 . 69444.1 

\ 

. 0 

0 

0 

0 J 





where x = 


6 

0 

6 

5 


and u 


T, 

T 2 




0.69444 0 

0 0.01 


(38) 


(39) 


We assumed perfect measurement and that plant disturbances enter the system coincident with the control 
torques. The sampling/update rates are given in Table 1 . 


TableJ^amgling/^xk^ 

Sample/Update Rate 

0 0.225 s 

6 0.028125 s 

T, 0.225 s 

0.028125 s 
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Five different compeasators were designed: an analog LQR, a multirate lead/lead, an optimized multirate 
lcad/lead, an optimized multirate general 2 nd order, and an optimized single rate general 2 nd order. We used a 
smooth step input to 5 re ^and 0^ defined as follows: 

fO.OOsf 1 - rn t<T c 

' c 

[0.01 m t>T c 

(40) 

Q rel it) = §rgy<r) , T c = 0.125 see 

p Li + L 2 

and the servo configuration shown in Figure 5 to measure the performance of the different compensators. The 
response of the TLA for the five compensators is shown in Figures 6a-6c. 

The analog LQR compensator used full state feed back. We provided this compensator as an example of 
the response possible using the cost function weighting matrices of (38). 

The multirate lead/lead was found using successive loop closures. We designed the control loops in the 
discrete domain so that the eigenvalues of the closed loop system matched those we obtained using LQR 
transformed to discrete time. This compensator consists of two simple lead loops: one from 5 to T 2 operating 
at the fast sampling/update rate, and one from 0 to Tj operating at the slow sampling/update rate. 

The final three compensators were synthesized using our new algorithm and the cost weighting matrices 
used to design the analog LQR compensator. The optimized multirate lead/lead was found by optimizing the 
pole/zero locations and gains of the lead/lead compensator found by successive loop closures. 

The optimized multirate general 2 nd order compensator uses the same sampling/update scheme as the 
lcad/lead compensators but has the compensator structure of (41), where a-, b tJ , c ijt and d- are the parameters 

which were optimized. This compensator has the maximum number of independent free parameters possible 
for a second order system |7]. 


a\\ 

0 

B = [ 1 fc 12 

c=\ Cn 

C\2 

D=\ d " 

d 12 

(41) 

0 

^22 . 

t>2l 1 J 

l C21 

Cll - 

- “21 

dn . 



The optimized single rate general 2 nd order compensator is a single rate equivalent of the multirate general 
2 nd order compensator. It has the same structure as the multirate general 2 nd order compensator, (41), but uses 
a single sampling rate. This sampling rate was chosen such that the number of computations required to 
implement either the multirate or single rate compensators during real-time operation are the same. 
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Our results are the same as those obtained using Yang’s algorithm. Thev demonstrate how multirale 
compensators can provide better performance than single rate compensators by trading lower bandwidth 
control of the slow modes for higher bandwidth control of the fast modes. In this example we were able to 
reduce the tip response over shoot 40% and the peak control torque 25% by using a multi rate controller over a 
single rate controller. 


Conclusions 

In this paper we have presented a new algorithm for synthesizing reduced order multirate compensators. It 
can be use to design compensators of arbitrary structure and dynamic order, with independent sampling/update 
rales for the compensator inputs, outputs and states. This algorithm provides the versitility of Yang’s aloruhm 
with out the numerical complexities associated with the finite time cost (unction. 

Finally, we do not want to discount Yang s algorithm altogether because, while our algorithm requires an 
initial stabilizing compensator, Yang’s does not. For those problems where finding an initial stabilizing 
compensator is difficult, we can always use Yang’s algorithm to find a stabilizing compensator and then 
quickly optimize the compensator parameter values with our algorithm. 
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APPENDIX 


Given a P k which stabilizes the multirate system we can calculate the steady state values of A k where A k is 
defined by equation (32) rewritten here as (A.1). 


0 = Qi + MP k H + (MP k Hf + (P k HfQ 2 P k H + F T ck A k + x F ek - A k 
for with A k = A k+N 


First simplify (A. 1) by defining 


' Qi 

M 

and F k = \ ^ 

. M T 

02 . 

IPkHl 


/ is an identity matrix 


(A.1) 


(A.2) 


Then (A.1) can be written as 


A* = F k Q$f k + F T ck A k+x F ck for k=l,2,... y N with A k = A k+N (A.3) 


Equation (A.3) represents a periodically time varying Lyapunov equation. We can create an equivalent 
single rate system by repeated application of (A.3). 

A* = JdkQdtdk + FhcA/cFdk for k= \,2,...JV with A k = A k+N (A.4) 


Fdk = Fc(k+N-\)Fc(k*N-2)Fc{k+N-3y • F ck 


(A.5) 


Fdk = 


F (*+N- \)Fc(k+N-2Fc{k+N-iy • F ck 
F(k+N- 2 )F c(k+N-3y ■ F ck 


Fk 


(A.6) 


Qd= 


03 0... 
0 Qy 


O' 

0 


0 0 003 
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Equation (A.4) is a time invariant Lyapunov equation which can be solved for A k . Once any A k has been 
found, the propagation equation (A3) can be used to find the remaining A k , 
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Figure 1. A General Multirate Compensator 



Figure 2. Example of a Multirate Sampling Scheme 
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Figure 3. Flow Chart of Optimization Algorithm 


78 












Figure 5. TLA Plant/Compensator Configuration 
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Figure 6a. Tip (8) Response to a Smooth Step Command to Tip Position 
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Figure 6b. Control Torque Tj Response to a Smooth Step Command to Tip Position 
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Figure 6c. Control Torque T 2 Response to a Smooth Step Command to Tip Position 
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I. Introduction 

There are many established methods for synthesizing multirate compensators [Berg, Yang, Mason, 
Glasson, Amit] but surprisingly few methods for analyzing the robustness of these systems. Current 
robustness analysis methods rely principally on the transfer function of the system. A multirate transfer 
function, in the traditional sense, does not exist, because most mulliratc systems are periodically time vary ing. 
Without some modification, established analysis methods cannot be applied directly to multirate systems. 

[Thompson] and [Apostolakis] have both proposed ways to extend existing robustness analys.s techniques 
to multirate systems. Thompson used “Kranc” operators to transform a special class of multirate systems, 
derived from sampled continuous systems, into MIMO single rate systems. Apostolakis transformed the 
general multirate system into a discrete time single rate MIMO system and then used impulse modulation to 
produce a continuous MIMO system [Boykin’s]. In both cases, the inputs and outputs of the new MIMO 
system were comprised of delayed samples of the inputs and outputs of the multirate system. Thompson and 
Apostolakis then used multivariable nyquist criterion and/or unstructured singular value analysis to calculate 
the gain and phase margins for the multirate system. 

In this paper we will present an alternative approach for extending nyquist criterion and singular value 
analysis to multirate and periodically time varying systems. Like Apostolakis. we transform the original 
system into an equivalent time invariant single rate system. However, we perform the robustness analysis in 
the “z” domain. By working in the “z” domain we can establish relationships between a mullirate/periodically 
time varying system and its time invariant single rate equivalent. These relationships clarify the limitations of 
nyquist and singular value analysis using single rate equivalent systems. 

The paper is divided into five section. Section I provides some back ground information about multirate 
systems and discusses transfer functions for multirate and periodically time varying systems. Section II 
discusses the application of the nyquist stability criterion to these systems; Section m discusses the 
application of structured and unstructured singular values analysis to these systems. Section IV contains a 

example of robustness analysis using structured singular values for a multirate system. Concluding remarks 
follow in Section V. 

I. Multirate and Periodically Time Varying Systems 

Before discussing robustness analysis we will first establish the relationship between multirate and 
periodically time varying systems. Then we will define an equivalent single rate system which will allow us to 

combine periodically time varying, mulnrate and time invariant systems using traditional block diagram 
techniques. 

A general multirate compensator is shown in Figure 1. Each input (y), output ( u ), and state (z) is 
sampled/updated at a rate which, in general, represents the desired bandwidth of the input or output with 
which it is associated, y is the value of y currently available to the digital processor from the zero order hold; 
while B is the current output from the digital processor which is held with a zero order hold to form the output 
u. A discussion of this compensator structure can be found in [Berg & Mason, and Yang]. 

Associated with this multirate compensator is a multirate sampling schedule which specifies the 
sampling/update rate for each input, output and state. We define the greatest common divisor of all the 
sampling/update periods as the shortest time period (STP) and the least common multiple of all the 
sampling/update periods as the basic lime period (BTP\ The integer N is defined as; 
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( 1 ) 


N = 


BTP 

STP 


When ihe sampling/update rates have ratios which are rational numbers, the sampling/update schedule is 
periodically time varying and the mullirate compensator can be modeled as a linear periodically time varying 
system of the form [Mason & Berg|: 

x(k+l) = A{k)x(k) + B(k)u(k) (2) 

yk) = C(k)x(k) + Dk)u(k) (3) 

where A(k) = A(k+N), B(k) = B(k+N), C(k ) = C(k+N ), and D(k) = D(k+N) 


The sampling period for (2)-(3) is one STP and the period of repetition is one BTP. 

Any practical multirate system can be modeled as linear periodically time varying system of the form of 
(2)-(3). Therefore, we will focus the remainder of the discussion on linear periodically time varying systems, 
of which multirate and single rate are a special case. 

Given a periodically time varying system of the form of (2)-(3). we can create an equivalent time invariant 
system by repeated application of (2)-(3) over the BTP [Meyer & Burrus]. The equivalent time invariant 
system is 


x(N(k+l)) = Ax(Nk) + Bu(Nk) 
y(Nk) = Cx(Nk) + Tfu(Nk) 


Where 


A = A(N-\)A(N-2)- ■ A(0) 


(4) 

(5) 

( 6 ) 


B = [A(N-\)A(N-2)-- A(\)B(0)\A(N-l)A(N-2) --A(2)B(l) | | B(N-l)] (7) 

C(0) 

C = C(1 ^ (0) ( 8 ) 

C(N-\ )A(N-2)---A(0) _ 


D = 


D( 0) 

C(1)B(0) 

C(2M(1)B(0) 

C(N-l)A(N-2)- ■ A(2)B(0) 
0 

D(l) 

C( 2)B(l)) 


C(N-\)A(N-2)- ■ -A(2)B(1) 


0 


0 


0 

D(N- 2) 

C(N-\)B(N-2) 


0 

0 

D(N- 1) 


(9) 


where y(Nk) = 

y(Nk) 

y(Nk+\) 

and u(Nk) = 

u(Nk) 
u(Ni t+1) 


j(Nk+N-\)_ 


u(Nk+N-\) 


86 



Equations (6)-(10) transforms the linear periodically time varying system, (2)-(3), with p inputs, q outputs 
and a sampling period of one STP to a linear time invariant system. (4 )-(5), with Np inputs, Nq outputs and a 
sampling period of one BTP. We will refer to (4)-(5) as the equivalent single rate system ( ESRS ) of (2)-(3). 
It is important to keep in mind that the inputs and outputs of an ESRS are comprised of samples of the inputs 
and outputs of a lime varying system. A consequence of this is the relationship: 

n 2 -bl 2 on 


We can calculate the transfer function of (4)-(5) using the following definitions for the Z Transform. Let 


ZUW)=£ x(i)z ‘ 


i~0 


Zyv(j:(A))=£ xiiN)!-™ 
1=0 


then 


y(z»,l) = Zsiyik+D) 

y(z N ) = 


r Mym 

ZaM*+ 1)} 


Z/v{y(*+AM)} 


with a similar definitions for u(z N ) and u(z N ,l). 
The transfer function for the ESRS, (4)-(5), is 


( 12 ) 

(13) 

(14) 

(15) 


y(z A 0 = G w (z A, )«(z A 0 

where G N (z N ) = C(fz N - A)' 1 /? + D 


(16) 

(17) 


The transfer funcuon is written as G N (z N ) to emphasize that the sampling period of the ESRS is one BTP , or 
N times the sampling period of the time varying system (2)-{3). 

So far, the ESRS has only been applied to periodically time varying systems. We could, however, 
calculate the ESRS of a time invariant system - in this case N can be any integer. 

The transfer function for the ESRS of a time invariant system can be calculate using (17). Alternatively, 
G n (z n ) can be calculated directly in terms of the transfer function of the time invariant system, G(z). Given 

y(z) - G(z)u(z) ' 18) 

and lollowing | Meyers & Bumis) we can write 


N - 1 

y(z)=£ r-'y(z*,/) (19) 

/=0 

AM 

y(r /v ,/) = ^-2 , X W) ( 2 °) 

N i = 0 
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Combine (18)-(20) to obtain 


?! 

where <p =e 14 


y(z N J) = ^rZ 1 ^ X 2 - m <p' ,{lm) G(z<l)‘)u(z N ,m ) (21) 

N 1=0 m = 0 

From (15) and (21) the I th row and rrf h column of Gn(z n ) is given by 

G N (z») l/n = ±z‘-”' X f^G^') (22) 

N hO 

For a time invariant system, G N (z N ) is made up of time and frequency shifted versions of G(z). A special 
case of (22) occurs when G(z<j>‘) = G(z ) for < = 0.1,.. ..AM. 

, |G(z) if l = m 

If G(z0 )=G(z) for i=0,l....A-l then G*(z "),,„ = | Q olherwise (23 > 

The simplest G(z) satisfying (23) is G(z) = constant. Equation (23) is an important relationship which will be 

used in Sections II and III. 

Equations (6)-(10), (16)-(17) or (21) can be used to compute state space and transfer function descnptions 
for the ESRS of a periodically time varying or time invariant system. The advantage of the ESRS is that is it 
allows us to manipulate time invariant and periodically time varying systems (e.g. multirate) as if they were 
both time invariant. The state space or transfer functions descriptions can be used to calculate input-output 
relations for systems in series or in a feedback loop just as in classical control [Khargonekar], In addition, 
[Kono] has shown that if the ESRS is stable then the time varying system from which it was derived will be 
stable. So, we need only worry about the stability of the ESRS. 

II. Nyquist Stability Criterion 

We can determine the stability of the periodically time varying system in Figure 2 by applying standard 
multiloop nyquist [McFarlane ...] criterion to the ESRS , since the periodically time varying system will be 
stable if its ESRS is stable. The ESRS return difference is given by 

/ - G^z^A^z") (24) 

and the nyquist contour is 

Z N = 0< w < 2 n 

When the periodically time varying system is SISO, we can determine traditional gain and phase margins 
from the nyquist plot. Recall that when A(z*‘) = A(z), A^z") is a diagonal transfer function matrix with A(z) 
on the diagonal. Thus, the ESRS for Figure 2 with gain and phase uncertainty can be written as 

G\{z^)actual ~ Gtf(z ^)nominatf NzN^^ ( 25 ) 

Phase and gain margins from the nyquist plot can be interpreted in the traditional sense even though the ESRS 
is MI MO because the inputs and outputs are correlated in time and a constant gain applies equally over all time. 
(Thompson] arrived at this same results using Kranc operators. 
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When ihe time varying system is MIMO. the standard MIMO nyquist restrictions apply. For MIMO lime 
varying systems, it is best to use a norm based approach such as singular value analysis. 


III. Structured and Unstructured Singular Values Analysis 

In the previous section we saw that the multiloop nyquist stability criterion can be applied to an ESRS to 
determine the stability of a periodically time varying system. In this section we will see that, with some 
limitations, both structured and unstructured singular value analysis can be applied to the ESRS to determine 
the robustness properties of a periodically time varying system. 

Given stable transfer functions G(z) and A(z) it has been shown that the closed loop system will remain 
stable through out continuous changes in A(z) if 

del (/ - G(z)A(z)) * 0 (26) 

or g (/ - G(z)A(z» > 0 (27) 

is satisfied around the nyquist contour, subject to certain restriction on Gz) and A(z) [Maciejousky,...]. By 
direct application of (26)-(27), a periodically lime varying system will be stable if 

del (l - G N (z")A N (z")) * 0 (28) 

or g(/ - G n (z n )A n (z n )) > 0 (29) 

is satisfied around the nyquist contour, because a periodically time varying system will be stable if its ESRS is 
stable. From (28)-(29) it follows that most singular value robustness tests can be applied directly to a ESRS to 
determine the robustness properties of a periodically time varying system. The results, though, must be 
interpreted in light of the fact that some of the inputs and outputs of the ESRS are time correlated. In the 
following paragraphs we will discuss the important differences between a single rate system and an ESRS and 
how these affect singular value analysis. 

1 ) There are qN, not q, singular values associated with each point on the nyquist contour for the ESRS of 
a time varying system with only q inputs and outputs. The additional singular values come from the time 
correlated inputs and outputs of the ESRS. Remember that the sampling period of the ESRS is one BTP, N 
times slower than the time varying system from which it was derived; but the ESRS has N times as many 
inputs and outputs as the original time varying system. The key point is that all of these singular values are 
important in determining the robustness of a periodically time varying system. 

If an ESRS is generated from a time invariant system, the singular values of the ESRS and the singular 
values of the original single rate system are related by the following expression. 

oGsieJ"*) = foG(0 V"). ctG(0 '«><“), . . . aG(<p" l en] (30) 

In (30), singular values associated with frequencies above l/BTP in G(z) are reflected back to lower 
frequencies in Gn(z n ). It follows directly from (30) that 
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||G / v(z A OlL=l|fi(z)IL 


( 31 ) 


where ||G(z)|L=sup olG(«>-)] 

2) The ESRS imposes a structure on any uncertainty A ^ Any A N in (28)-(29) must obey (17) or (22); 
this automatically imposes a structure on Ayy.. The problem is, we are often interested in a time invariant plant 
uncertainty. A, that destabilizes the system and not in A N . a(A w ) found using unstructured singular value 
analysis is often overly conservative because it accounts for not only the fictitious perturbations normally 
associated unstructured singular values but also for lime varying and non-causal perturbations. It is important 
to remember that C(A N ) found using unstructured singular values can be extremely conservative and may not 
reflex (7(A). 

Structured singular values provides a mechanism for finding A. For an ESRS with a time invariant 
uncertainty, A(z), the definition of the structured singular value, ft. can be written as 

JO if det(I - G\(z n )Ah(z n )) * 0 for any A e A B d ^ 

V- (G N (zN )) = |( min (o(A(z)): del (/ - G N (z)A N (z)) = o) ' otherwise 

A bd is the form of the permissible block diagonal perturbations A; and the structure of A N must satisfy 
equation (22). 

Unfortunately for a general A(z), the structure of A^z) is often very complex and finding a good estimate 
of size of A is difficult. However, when A is a constant, as is the case for many problems, 

A jy = diag( A, A, ... A) with N blocks (33) 

and structured singular value analysis can be used to determine A. 

When A is a time varying, but not a function of “z”, A N becomes 

A n = diag( A(l), A(2), ... A (N)) (34) 

and has no repeated blocks. Equation (34) must be interpreted with care - (34) implies that the value of A(k) is 
constant over the sampling interval, STP, and changes instantaneously to A(*+l) at the next sampling instant. 
This may not be a good model of time varying uncertainty. 

3) When A is a constant then each A block of A N can be scaled independently. Using the block diagonal 
scaling properly of |i (Maciejowski|, and (33) or (34) it is straightforward to sec that 

If A is a p by p matrix then ^(DGsit^O '*) = p^G/v^z^O) (35) 

where D = ( d x I p , ... Vp) and / is a p by p identity matrix 

In addition if A is block diagonal then each of sub-block of A can be scaled in a similar manner. 

An interesting result of (35) is that the upper bound for p(G^/(z^ 0) given by 
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\l(G n {z n )) < itfiDGsiz^D 1 ) 


(36) 


is the same whether A is periodically time varying or time invariant, subject to the interpretation of a time 

varying A mentioned in item 2. The upper bound for the time invariant case found using (36) is of course 
more conservative. 


4) Singular value plots of ESRS transfer function matrices should not be interpreted in the frequency 
domain. The ESRS has time correlated inputs and outputs. The response from one input to one output 
represents only part of the total signal between the input and output of the periodically lime varying system 
A meaningful quantity for an ESRS is its infinity norm. From (1 1) and [Francis] we can write that 

lt>l 2 IMI 2 „ .. 

SUP \W =SUP W = ,|Gw(Z ^ for l|u|l 2 (37) 

2 2 

Thus, HGMz^L can be interpreted as the maximum gain of the system for all u with a bounded two norm, 

just as in the single rate case. For the single rate case the maximum gain occurs when u is sinusoidal - this is 
not necessarily true for the periodically time varying system. 


Singular value analysis of an ESRS, both structured and unstructured, can be used to determine the 
robustness of periodically time varying system. As we have discussed, there are limitations to this analysis 
because the inputs and outputs of an ESRS are time correlated. 

IV. Two Link arm Example 

The results of the previous sections will be illustrated by calculating the gain margins for a planar two link 
robot arm (TLA) using structured singular values. Two different cases are considered: 1) the TLA with a 2 nd 
order multirate compensator and 2) the TLA with a 2^ order single rate compensator. 

The TLA is shown in Figure 3 and is described further in [Berg and Yang], The two compensators were 
designed to minimize a cost function quadratic in the states and controls using the optimization method 
described in [Mason & Berg]. 

The 2 nd order multirate compensator uses the compensator structure of (38), where a r , b r , c r , and d t are 
the parameters which were optimized. 


A = 


ait 

0 


0 

322 



b)2 

1 


c = 

cm 

C ' 2 1 D = 

d„ 

d|2 


L C21 

C 2 2 J' 

d2i 

d22 . 


(38) 


The sampling/update rates for the compensator are listed in Table 1. In addition the compensator state 
associated with 0 and T, is updated at the slow rate while the state associated with 8 and T 2 is updated at ihe 
fast rate. For the multirate compensator STP = .028125, BTP = .225 and N = 8 
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““ 

Mullirate Compensator 

Single Rale Compensator 

0 

0.225 s 

0.05 s 

8 

0.028125 s 

0.05 s 

T l 

0.225 s 

0.05 s 

t 2 

0.028125 s 

0.05 s 


Table 1. Sampling/Update Periods for the Compensators 


The single rate 2 nd order compensator is the single rate equivalent of the mullirate compensator. It has the 
same structure as the 2 nd order multirate compensator and minimizes the same cost function, but uses a single 
sampling/update rate. This sampling/update rate was chosen such that the number of computations required to 
implement either the multirate or single rate compensators during real-time operation is the same. The 

sampling/update rate for the compensator is shown in Table 1 . 

A block diagram of the TLA, compensator, and output gain uncertainty is shown in Figure 4. The bock 
diagram in Figure 4 can be cast into the standard structured uncertainty model shown in Figure 5 where the 
gains k x and ^ are allowed to vary independently. An upper bound on the structured singular values for the 

multirate and the single rate cases was calculated using the following [Safonov], 

p(Q) < ir^a(DQP~ x ) ^ X./>(G) 
where \ P (Q) is the Perron-Frobenius eigenvalue of Q 

For the multirate case, the ESRS for the plant, compensator and uncertainty was calculated for A = 8. 
They were combined as shown in Figure 5 and an upper bound for p. as z traversed the nyquist contour, was 
calculated using (39). A lower bound on the maximum singular value of the gain matrix is given in Table 2. 

For the single rate case an upper bound on p was calculated using two different methods. First, p was 
calculated directly for the single rate system, N = 1 . An exact value for p can be calculated because there are 
only two blocks in the uncertainty matrix [Doyle]. Next, an ESRS was constructed for the single rate case 
using N = 8. An upper bound for p was calculated using (39). For the ESRS system, the uncertainty matrix 

has 8, 2 by 2 blocks. These results are summarized in Table 2. 


Design 

Gain Margin 0 

- — 

*1 oi 
_0t 2 . 

_ i_ 

Multirate 2 nd Order 

Single Rate 2 nd Order 

Single Rate 2 nd Order using ESRS 

0.535 

0.513 

0.389 


Table 2. Gain Margins for TLA with Multirate and Single Rate Compensator 


The two p estimates for the single rate case illustrate the disadvantage of using (39) to calculate the upper 
bound of p for an ESRS systems. As in the multirate case, the ESRS single rate case accounts for periodically 
time varying uncertainties, resulting in a conservative estimate of p. See item 3, Section III. 
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For the assumed uncertainty model the multirate compensator was slightly more robust, even given the 
conscrvativcness of the estimate for p. The multirate compensator is able to compensate for larger gain 
uncertainty because it has higher bandwidth control of the second link than does the single rate compensator. 


V, Summary and Conclusions 

In this paper we have shown how nyquist criterion and singular value analysis can be applied to multirate 
and periodically time varying systems using their ESRS. For SISO systems, traditional gam and phase 
margins can be found by direct application of the nyquist criterion to the ESRS. For MIMO systems, 
structured singular values can be used to determine the maximum size of an uncertainty. The results of 
singular value analysis, though, must be interpreted in light of the fact that some of the inputs and outputs of 
the ESRS are time correlated. We pointed out several important resulting limitations of singular value analysis 
using an ESRS. Finally we demonstrated robustness analysis for a two link arm with a multirate compensator 
using structured singular value. 
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APPENDIX 

Lemma: oGv(e> N<0 ) = [oGty'W oG(<p V"). • • • oG(<p N 'e'®)] 
Proof: From the definition of a transfer function we can write 

y(z) = G(z)u(z ) 



y(0°z) 


«(0°z) 

where y(z) = 

y(<p x z) 

, «(z) = 

«(</>' z) 


L >'(0 z )_ 


,,N-1 ^ 
_U(0 z)J 


G(z) = 


G(<p z) 0 
0 G(<p *z) 


0 0 G(0 W 'z) j 


(A.l) 


(A. 2) 


AM » 

From (12), >(z) = £ z-^(z N ,/) so that we can write y(z) = Ty(z ) 
/=o 


where T = 


/ z-'/ 


Z -(JV-1)/ 

p'z)< N -"l 


(A. 3) 


7 has the property that 7T* = 2V/ if z is evaluated on the unit circle and / is an identity matrix of 
appropriate dimensions. 

Then _ 

y(z*) = r'G(z)Tu{z N ) so that G N (z N ) = T x G{z)T (A.4) 

Now using the fact that a 2 (A) equals the eigenvalues of A*A. and that the eigenvalues of a block diagonal 
matrix arc the eigenvalues of the individual block it follows that 


oGn^") = [cG(0 V"), aG(<p X e> a ), ■ • • oG(<p N l ei b> )\ (A.5) 
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Figure 1 . A General Multirate Compensator 



Figure 2. Periodically Time Varying System 
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l 2 



Parameters: Mass Length 

L , 1.235 kg 0.965 m 

L 2 0.163 kg 0.167 m 

Inputs: Torque Tj and T 2 
Outputs: 9 and 8 


Figure 3: Diagram of the Two Link Arm 
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Figure 4. TLA Feedback Loop with Output Gain Uncertainty 



Figure 5. Structured Uncertainty Model for TLA 
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